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Abstract 

The  report  presents  heat  transfer  modeling  of  Thrust  Vector  control 
systems  using  the  PHOENICS  computer  code. 

Simple  two-dimensional  wedge  and  blunt  bodies  have  been  examined  in 
supersonic  cold  flow,  for  both  laminar  and  turbulent  flow  cases. 

The  research  presents  a  numerical  solution  of  the  supersonic  compressible 
viscous  two-dimensional  flow  field.  Post  calculations  were  done  to  estimate 
skin  friction  coefficient,  surface  heat  flux,  heat  transfer  coefficient  and 
Stanton  number  distributions  in  both  wedge  and  blunt  cases. 
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1 .  Introduction 


This  report  describes  a  numerical  analysis  of  heat  transfer  of  a  typical 
jet  vane  configuration  used  for  thrust  vector  control.  The  work  was  carried 
out  under  contract  Nos.  N6227 1-85-M-0443  and  N6227 1-86-M-0206 ,  for  the  Naval 
Postgraduate  School. 

The  tasks  to  be  accomplished  under  the  first  contract  were: 

Task  I:  Formulate  the  conservation  equations  of  momentum  energy  for 
two-dimensional,  supersonic  flow  in  geometries  typical  of  thrust  vector 
control  systems. 

Task  II:  Formulate  boundary  conditions  for  these  equations  appropriate 
to  thrust  vector  control  systems. 

The  tasks  to  be  accomplished  under  the  second  contract  were: 

Task  I:  Continue  and  update  the  formulation  of  thrust  vector  control 

geometries  based  on  the  input  from  the  Naval  Weapons  Center 
( NWC) • 

Task  II:  Construct  the  computational  model  for  implementation  in  the 

PHOENICS  code,  of  the  thrust  vector  control  geometries  and 
flow  conditions  provided  by  NWC. 

Task  III:  Run  the  PHOENICS  code  for  the  previously  formulated  models. 

Analyze  and  interpret  the  PHOENICS  results  for  surface 
temperature  *and  heat  flux. 

Thrust  vector  control  components  such  as  jet  vanes  and  jet  tabs  are 
exposed  to  high  speed  hot  gases  at  the  exit  of  a  rocket  nozzle. 

Estimation  of  the  heat  transfer  from  the  hot  exhaust  gases  to  the  vane  is 
major  consideration  in  the  correct  design  of  a  vane,  and  its  ability  to 
survive  during  its  mission. 
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The  research  work  was  done  under  the  framework  of  the  tasks.  A  brief 
survey  of  what  has  been  done  according  to  the  task  is  given: 

Task  1  (M-0443):  Heat  transfer  modeling  of  thrust  vector  control  vane 
requires  supersonic  compressible  viscous  flow  analysis. 

In  order  to  meet  the  requirements,  the  conservation  differential 
equations  of  mass  momentum  energy  and  the  two  k-e  turbulent  equations  were 
formulated,  and  additional  algebraic  formulas  for  the  relations  between 
pressure  density  and  the  equation  of  state  for  ideal  gas. 

Task  II  (M-0443):  The  physical  dimensions  of  the  flow  field  grid  were 
chosen  and  the  boundary  conditions  for  the  Navier-Stokes ,  energy  and  the  two 
k-q  turbulent  model  equations  were  given. 

Task  I  (M-0206):  Working  on  the  task,  the  actual  configuration  of  a  jet 
vane  that  is  presently  being  tested  at  NWC  has  been  modeled.  The  geometry 
being  used  is  a  wedge  which  has  the  same  half  angle  and  dimensions  as  the  NWC 
jet  vane. 

Task  II  (M-0206):  BFC  (Body  fitted  coordinate)  version  of  PHOENICS  code 
(Ref.  3)  was  used  for  calculating  the  flow-field  and  heat  transfer  over  the 
model.  Using  the  BFC,  a  better  geometrical  approximation  to  vane  shape  could 
be  achieved. 

Non-Uniform  grids  have  been  utilized  in  order  to  model  complicated  regions 
in  the  flow  field.  Relaxation  parameters  and  false  times teps  options  were 
adjusted  to  enable  efficient  computer  runs  with  good  convergence. 

Task  III  (M-0206):  In  carrying  out  this  task,  four  major  runs  have  been 
analyzed : 

Two  geometric  configurations  were  used:  wedge  vane  and  blunt  vane  (see 
Figures  1,  2,  3,  4);  each  one  in  both  laminar  and  turbulent  flow  conditions. 
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Numerical  results  for  fluid  field  and  thermodynamic  properties  of 
pressure,  temperature,  density,  Mach  number  and  velocities  are  given  in 
appendix  C. 

Post-calculations  of  heat  transfer  coefficient-,  skin  friction  coefficient 
and  Stanton  number  are  given  in  Figures  (6,  7,  8,  9,  10,  11). 

The  next  chapters  describe  in  more  detail  the  process  of  building  the 
model  and  the  analysis  of  the  results. 


2.  PHQENICS  Description 


The  present  work  addresses  the  heat  transfer  modeling  of  thrust  vector 
control  systems*  In  this  effort  the  Navier-Stokes  approach  is  applied  by 
using  a  computer  code  which  is  capable  of  simulating  a  large  number  of  fluid 
flow,  heat  transfer  and  chemical  reaction  processes  which  arise  in  industry 
and  elsewhere.  This  code  is  called  PHOENICS,  which  is  an  acronym  standing 
for:  'Parabolic,  Hyperbolic  or  Elliptic  Numerical  Integration  Code  Series.1 

The  name  comes  from  the  fact  that  the  differential  equations  of  fluid  flow, 
etc.  arise  in  forms  classified  by  mathematicians  as  parabolic,  hyperbolic  or 
elliptic;  and  PHOENICS  solves  these  equations,  whatever  their  form. 

Built  into  PHOENICS  are  the  major  conservation  laws  of  physics  (mass, 
momentum,  and  energy)  applied  to  a  large  number  of  continuous  subdomains 
called  'cells,1  into  which  the  domain  of  study  is  artificially  divided.  The 
number  of  cells  can  be  few  or  many  according  to  the  requirements  of  the 
problem.  Because  of  numerical  stability  considerations  the  restrictions  on 
cell  refinement  can  become  particularly  burdensome  in  the  calculation  of  a 
turbulent  boundary  layer  where  a  very  fine  mesh  near  the  wall  may  be 
required. 

When  supplied  with  appropriate  information  concerning:  the  physical 
properties  of  the  materials,  the  geometrical  and  other  constraints,  the  inlet 
and/or  initial  conditions,  PHOENICS  computes  the  corresponding  solutions  to 
the  relevant  differential  equations,  expressing  them  as  tables  of  numbers 
describing  the  field  of  velocity,  temperature  concentration,  etc. 

Detailed  information  about  PHOENICS  is  given  in  [Ref.  3]. 
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2.1  The  Structural  Principle  of  PHQENICS 


The  code  consists  of  three  major  parts:  Satellite  subroutine, 
Ground  subroutine  and  Earth  library. 

The  satellite  subroutine  is  the  main  input  subroutine  and  should 
provide  the  answers  to  the  questions: 

-  what  kind  of  process  is  to  be  simulated 

-  what  are  the  properties  of  the  fluid 

-  what  are  the  shape  and  size  of  the  domain 

-  how  fine  is  the  grid  to  be  employed 

-  to  what  degree  of  accuracy  is  the  calculation  to  be  continued 

-  and  what  output  should  be  provided 

Ground  subroutine  is  active  during  the  computing  process  and  is  used  for 
updating  properties  which  vary  with  time,  temperature,  etc.  For  example: 
viscosity  depends  on  temperature  or  density  depends  upon  pressure  and 
temperatures ,  etc • 

Earth  library  is  the  main  solver  generator.  It  is  given  as  a  binary 
library  and  does  not  enable  the  user  access  to  the  source  code. 

2 . 2  Numerical  Scheme 

The  numerical  scheme  used  by  the  code  is  the  simpler  (semi- implicit 
method  for  pressure-linked  equations  revised)  (Ref.  9).  The  scheme  was 
developed  by  Patankar,  S.  V.  and  Spalding,  D.  B. 

The  scheme  requires  an  additional  dependent  variable,  the  pressure 
correction,  which  has  no  physical  meaning  but  should  take  part  in  the 
process . 

The  value  of  the  pressure  correction  should  tend  to  zero  in  the 
convergence  process • 

Two  additional  differential  equations  are  solved:  for  the  pressure,  and 
for  the  pressure  correction. 
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3 .  Geometry  and  Dimensions 


Symmetrical  2-D  planar  geometry,  which  is  shown  in  Figure  2,  was  chosen 
to  be  the  approximation  of  the  MWC  vane  in  Figure  1. 

Two  geometrical  profiles  were  examined,  one  with  wedge  leading  edge  and 
the  second  with  blunt  leading  edge. 

The  dimensions  of  the  domain  in  Figure  3  and  4  satisfy  aspect  ratio  of 
10:1  in  the  vertical  y  coordinate,  A  high  aspect  ratio  in  the  coordinate  is 
important  for  the  assumption  of  free  stream  conditions  at  the  upper  boundary. 
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4.  Assumptions 


Postulating  the  right  or  the  wrong  assumptions  has  the  most  influence  on 
modeling  process.  The  stage  was  carried  out  very  carefully  in  order  to  make 
the  most  compatible  model  with  reality. 

4.1  Steady  state: 

The  modeling  assumes  steady  state  physical  phenomenon  process. 

— —  (all  properties)  *  0 
0  t 

This  is  a  valid  assumption  since  the  time  constant  for  the 
convection  process  is  much  shorter  than  the  time  constant  for  the  wall 
conduction. 

By  assuming  the  wall  temperature  to  be  constant,  the  two  procedures 
are  decoupled. 

In  hot  flow  it  is  important  to  run  the  code  for  a  wide  range  of  wall 
temperature  which  will  take  into  account  the  influence  of  different 
temperatures  on  the  heat  convection  process. 

4.2  Cold  Air  Flow 

Ambient  temperature  air  flow  which  was  utilized  by  NWC  experiments 
is  being  used  in  the  computations. 

4.3  Ideal  Gas 

The  gas  is  assumed  to  satisfy  the  ideal  gas  equation  of  state 

p  *  p  RT  (4.1) 

This  is  a  fairly  good  assumption  for  nonreactive  gas  flow.  In  spite  of  the 
values  of  static  temperature  can  decrease  to  200[k],  the  density  remain 
relatively  low. 
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This  assumption  is  an  important  simplification  to  the  solution  in  Ref.  10 
which  used  the  isentropic  relation  between  pressure  and  density  instead 


4.4  Constant  Pr,y: 

Prantdl  number  and  y  (ratio  of  specific  heats)  were  found  to  have 
negligible  variations  in  the  temperature  range  of  the  model.  (200k  t  350k) 

4.5  Varying  Viscosity  and  Thermal  Conductivity: 

Vi  and  k  are  much  more  dependent  on  temperature  especially  very  close  to 
the  solid  wall  where  values  of  y  and  k  influence  strongly  the  shear  and  heat 
transfer  mechanism.  To  account  for  the  temperature  dependence  power  law 
relations  have  been  formulated  for  y  and  k. 

»r  0.666 

u  ■  u0(f — )  (4.3) 

Ao 

7  0.666 

k  -  k  — )  (4.4) 

o 

4.6  Parallel  Flow 

Gas  flow  at  the  exit  of  the  exhaust  nozzle  is  more  likely  to  be  a  conic 
source  flow  than  parallel  flow. 

If  the  half  angle  of  the  nozzle  is  small,  (ot  <  15°),  parallel  flow  is  a 
good  assumption 

4.7  Negligible  Radiation 

Assessments  that  were  done  showed  that  heat  convection  is  at  least  one 
order  of  magnitude  greater  than  heat  flux  by  radiation. 
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4.8  Laminar  and  Turbulent  Solutions 


In  order  to  overcome  lack  of  ability  to  predict  transition,  separated 
laminar  and  turbulent  calculations  were  done  for  each  case.  The  turbulent 
solution  utilizes  the  (k-e)  eddy  viscosity  model  Ref.  5. 

4.9  Constant  Wall  Temperature 

The  vane  wall  is  assumed  to  have  constant  temperature  during  the  time  of 
calculation. 
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5 •  Governing  Equations 


The  conservation  equations  for  the  compressible  flow  of  the  mathematical 
model  consists  of  a  viscous,  Newtonian  perfect  gas  consisting  of  the  following 
six  differential  equations: 

Conservation  of  Mass: 


p)  +  v*(p$) 


0 


(5.1) 


Conservation  of  momentum: 


Cp  4>)  +  V*  -  uv<t>)  vp 


(5.2) 


where  <}>  is  V  or  W  velocity  component  for  y  and  z  direction. 

Conservation  of  Energy 

(ph)  +  V*(ptfh  -±-  Vh)  -•§£  (5.3) 

where  h  is  the  total  enthalpy, 
h  -  Cp  T0 

where  To  is  the  total  temperature 
To  -  Tstat*(l  +  ^M2) 


In  the  case  of  laminar  flow  the  governing  equations  (5.1),  (5.2),  (5.3) 
are  sufficient  to  determine  a  solution  when  proper  boundary  conditions  are 
applied  and  the  equation  of  state  (4.1)  is  provided. 

Turbulence  Model: 

In  turbulent  flow  it  is  necessary  to  hypothesize  a  turbulence  odel 
relating  the  turbulent  viscosity  to  the  other  problem  variables. 
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The  model  used  in  PHOENICS  is  the  eddy  viscosity  (k-e)  model  (Ref.  3, 
Ref.  5J.  In  this  model  k,  the  turbulent  kinetic  energy  and  e,  the  turbulence 


dissipation  rate,  are  treated  as  properties  of  the  flow  and  conservation 
equations  are  postulated  for  these  properties.  The  two  conservation  equations 
are:  one  for  k,  the  kinetic  energy  of  turbulence: 


Second  equation  for  e,  the  dissipation  rate  of  turbulence 


De  m  9  ,  veff  9 e 
Dt  *  9Xj '  oe  3Xj 


H-)  +|(ClGk-  C2e) 


(5.5) 


where 


3Ui  3  Ui  3Ui 


(5.6) 


(5.7) 


cj,  C2 ,  a^,  are  empirical  constants  which  are  provided  in  PHOENICS. 

The  reason  for  using  the  (k-e)  model  is  because  it  is  the  most  verified 
model  for  engineering  applications.  It  combines  simplicity,  universality,  and 
realism  of  predictions  in  most  cases. 

Two  additional  differential  equations  are  solved  also  in  order  to  satisfy 
the  SIMPLER  algorithm  as  was  mentioned  in  chapter  2.2.  The  description  of  the 
pressure  and  pressure  correction  equations  is  provided  by  Ref.  9. 
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6.  Input  Variables 


The  properties  of  mach  no.  stagnation  presence  and  temperature  of  the  gas 
were  provided  by  NWC;  additional  properties  were  taken  from  air  tables: 

Mach  number: 

Moo  ■  3*2 

Stagnation  pressure: 

Po  =  55. 105  [Pa] 

Stagnation  temperature: 

T0  =  555.55  [K] 

Gas  constant 
R  =  287.  [ J/kg*k] 

Specific  heat  ratio 
y  -  1.35 

Laminar  Prandtl  Number 
Pr  -  0.7 

Turbulent  Prandtl  Number 
Prt  -  0.9 

Constant  Pressure  Specific  Heat 
Cp  -  R/(1-1/y)  U/kg-k] 

Laminar  Viscosity 
V  -  0.1716  *  10-5  *  (T/273.)0*666 
Thermal  Conductivity 
k  -  p  Cp/Pr 

The  gas  properties  in  the  inlet  boundary  are  equivalent  to  the  properties 
at  nozzle  exit.  Inlet  properties  are  calculated  from  the  stagnation  values  in 
the  combustion  chamber.  The  calculation  was  done  by  assuming  one  dimensional 
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isentropic  expansion  from  combustion  chamber  to  the  nozzle  exit  (inlet  for 
the  vane). 


Pressure 

Pi  =  PQ/(  1  +  JCl  M2)Y/Y-1 

(6.1) 

Temperature 

Ti  =  T0/(l  +  M2) 

(6.2) 

Density 

Pi  =  P/RT 

(6.3) 

Enthalpy 

hi  =  STi 

(6.4) 

Sonic  Velocity 

Ci  =  /Wl 

(6.5) 

Velocity 

Wi  =  C*M 

(6.6) 

The  subscript  i  signifies  inlet  property 
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7 .  Boundary  Conditions 


The  flow  field  described  in  Figures  3,  4  has  four  boundaries,  which  can 
be  named:  inlet,  outlet,  frees tream  boundary  and  solid  wall. 

Super  sonic  flows  have  a  hyperbolic  mathematical  nature.  The  field 
consists  of  influence  zones,  the  flow  at  every  point  is  governed  only  by  its 
influence  zone,  basically  by  the  upwind  stream. 

As  a  consequence  from  the  discussion,  it’s  obvious  that  the  outlet 
boundary  condition  has  no  influence  on  the  upstream  flow.  The  boundary  values 
that  are  given  at  the  outlet  are  to  satisfy  some  numerical  needs  only. 

7.1  Inlet 

Parallel  uniform  flow  with  known  velocity,  enthalpy,  pressure,  and 
density:  equation  (6.1),  (6.3),  (6.4),  (6.6)  are  given  at  the  left  boundary 
of  the  grid.  In  PHOENICS  this  is  specified  as  the  LOW  side  of  the  first  Z 
cell  • 

In  turbulent  flow,  boundary  conditions  are  supplied  for  k  and  e.  The 
values  that  are  given  are  based  on  empirical  values: 

k±  -  0.0  W*  (7.1) 

e±  «  0.16  kl*S/(5*GH)  (7.2) 

where  GH  is  half  the  vane  thickness 

7.2  Outlet 

As  was  mentioned  previously,  the  outlet  has  negligible  effect  on  the 
results.  The  only  property  that  is  specified  at  the  outlet  is  the  pressure. 
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7 .3  Frees tream  Boundary 


Assuming  that  the  upper  boundary  is  chosen  to  be  far  enough  away,  the 
default  boundary  condition  option  of  PHOENICS  is  used.  This  implies  a  line  of 
symmetry  where  all  gradients  are  zero, 

7.4  Solid  Wall 

Zero  velocity  and  constant  wall  enthalpy  (temperature)  are  assumed  on  the 
wall.  In  PHOENICS  the  wall  is  the  SOUTH  side  of  the  first  y  cell.  The  high 
enthalpy  and  velocity  gradients  near  the  wall  demands  a  refined  grid  close  to 
the  wall.  Values  of  shear  stress  and  heat  flux  are  calculated  to  first  order 
accuracy  using: 


aw  Wi 

Tw  -  u  —  =  4 


(7.1) 


qw 


-H_  ill  ~  4_  •  hi 

Pr  3y  Pr  Ayi/2 


(7.2) 


In  turbulent  flow,  a  wall  function  is  used  to  provide  the  wall  condition 
for  velocity,  enthalpy,  k,  and  e 


7 .3  Wall  Function 

The  wall  problem  in  the  numerical  computation  of  flows,  especially  in 
turbulent  flow,  is  an  old  one  and  most  authors  have  adopted  similar 
techniques.  In  effect  they  "bridge  over"  the  region  very  close  to  the  wall  by 
introducing  special  functions  which  are  called  wall  functions.  These  are 
often  empirical  in  origin.  Accounts  may  be  found  in  Ref.  11. 

The  problem  arises  as  follows.  Turbulence  dies  out,  close  to  the  wall, 
because  the  no  slip  condition  and  the  rigidity  of  the  wall  make  all  the 
velocity  components  fall  to  zero.  The  consequence  is  that  the  effective 
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viscosity  and  other  transport  properties  fall  there  to  their  laminar  values 
and  the  result  is  a  rapid  variation  with  distance  from  the  wall  both  of  the 
$  '  s  and  of  their  gradients • 

Where  <j>  signifies  general  dependent  variable ,  It  is  possible  to  compute 
these  variations  in  detail,  by  using  a  computer  code  such  as  PHOENICS  on  two 
conditions: 

(i)  the  grid  points  must  be  packed  into  the  region  of  steep  gradient 
changes  closely  enough  for  sufficient  numerical  accuracy  to  be  obtained 

(ii)  the  functions  appearing  in  the  turbulence  model  equations  must 
properly  represent  the  influence  of  local  Reynolds  number  on  turbulence. 

Under  the  conditions  above,  the  wall  function  sequences  in  the  program 
act  as  follows: 

The  Reynolds  number  is  first  evaluated,  based  on  the  resultant  velocity 
parallel  to  the  wall,  on  the  distance  from  the  wall  to  the  grid  node  and  on 
density  and  laminar  viscosity.  If  this  Reynolds  number  is  less  than  132.25 
(the  value  at  which  the  laminar  and  turbulent  wall  function  intersect)  a 
laminar  wall  function  is  used.  If  this  Reynolds  number  turns  out  to  be 
greater  than  132.25  the  velocity  variation  is  logarithmic  and  the 
corresponding  shear  stress  coefficient  is  evaluated.  This  corresponds  to  the 
commonly  used  "log  law"  wall  function.  [Ref.  4] 

7 .6  Boundary  Conditions  in  Phoenics 

PHOENICS  utilizes  source  terms  for  creating  boundary  conditions.  The 
form  of  the  source  term  of  each  dependent  variable  <|>  is: 

S -  ( [m]  +  C*)  (V+  -  *p)  (7.3) 

where:  m  -  is  mass  flux  source 

<{>p  -  is  the  value  of  the  dependent  variable  at  point  near  the  boundary 

C<j>,  V -  two  coefficients  specified  by  the  user.  The  source  term  for 
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mass  flux  is  simply 

%  -  C.  <v«  -  V  <7-4> 

where:  Pp  -is  the  pressure  near  the  boundary  and  C*.  Vm  are  two  coefficients 

The  values  of  C*  and  V*  for  the  dependent  variables  in  SATELLITE  are:  At  the 


Inlet : 

C  =  2  T  — 

Lm  4  y=T  Wi 

Vm  “  PD  Pi/Po 

Cw  =  ch  =  Ck  =  Ce  =  °* 

Vw  =  w± 

Vh  “  hi 

Vk  =  Kj. 

VE  -  si 
At  the  Outlet: 

Cm  -  1000  *  Wi-pi/Pi 

Vm  ’  Pi 

At  the  Wall  (laminar) 

Cw  =  u/(0.5Awi) 

Vw  “  ° 

Ch  =  u/Pr/(0.5*Ayi) 

Vh  =  Cp*Tw 

At  the  Wall  (turbulent) 

Cw  ‘  Ch  '  ck  -  CE  -  WALL 

V„  -  Vk  -  V£  -  0 
Vh  =  Cp*Tw 


(7.5) 

(7.6) 

(7.7) 

(7.8) 

(7.9) 

(7.10) 

(7.11) 

(7.12) 

(7.13) 

(7.14) 
(7. is) 

(7.16) 

(7.17) 

(7.18) 

(7.19) 

(7.20) 
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8.  Mesh  Generation 


In  this  work  a  two-dimensional  mesh  is  being  used  with  18  x  29  cells  in 
the  y  and  z  coordinate  respectively.  A  Nonuniform  grid  has  been  used  for 
both  directions.  Figures  3  and  4  shows  the  grid  in  the  z  direction.  A  finer 

grid  is  used  in  the  blunt  region,  IZ  =  (7  *  17),  and  in  the  zone,  where  the 

inclined  wall  transitions  to  a  straight  wall,  IZ  -  (23  *  26). 

In  the  y  coordinate,  except  in  the  boundary  layer  region,  the  grid  is 

uniform.  To  obtain  a  finer  grid  resolution  in  the  boundary  layer  for  the 
laminar  flow  case  the  first  five  cells  in  the  y  direction  from  the  wall  obey 
the  following  proportionality  relationship: 

tv  3  Anax 

BY F RAC  (IY)  -  (-y-)  (8.1) 

Where  BYFRAC( IY)  is  the  distance  from  the  south  side  to  the  north  side  of  the 
cell  of  particular  interest,  divided  by  total  length  of  the  domain,  IY  is  the 
cell  number,  Amax  *s  roaximum  allowable  cell  height,  and  GH  is  the  half 
thickness  of  the  TVC  jet  vane. 

A  fine  grid  resolution  for  the  turbulent  flow  case  is  set  up  in  the  same 
way  as  laminar  flow.  The  only  difference  comes  from  the  selection  of  the 
first  five  cells  in  y  direction.  The  following  calculation  shows  the 
difference • 

From  the  laminar  solution  and  the  given  properties  the  following  are 
known: 

w  -  885.2[m/s] 

Ulam  “  1*10"5  [N.s/m] 

Po  -  5.5*106[Pa] 

Pstatic  -  1.048*10 5 [ Pa] 

Y  -  1.35 
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p  =  1.835  [kg/m] 

Using  the  values  above  and  the  length  of  vane,  which  is  0.095m,  A 
corresponding  Reynolds  number  was  calculated: 

Rej  .  i  -  <1^888^09S)  ,  ,.54  .  106 

Using  a  power  law  correlation  for  the  boundary  layer  thickness: 

j  =  0.37  *  Rez_l/5  (8.2) 

From  equation  (8,2)  the  boundary  layer  thickness  at  the  high  end  of  the 

domain  has  been  calculated  as  6  «  2  *  10“3[m] 

With  Re  based  on  WM  the  velocity  parallel  to  the  wall,  —2*  the  distance 


from  the  wall  to  the  first  grid  node,  the  density,  and  \*iam  the  laminar 
viscosity,  Ay  must  satisfy  the  condition 


^eA 


P»  Woo  Ay 


>  132.25 


or 


2  Ulam 

Therefore  the  interval  of  Ay  is  chosen  such  that 


AY  >  6.48  *  10“6 [m] 


2  *  l(T3[m]  >  AY  >  6.48  *  10“6[m] 


In  this  effort  using  the  relationship 
BYFRACt  IY)  -  (-S-)2  <£=§) 

Ay  has  been  calculated  as  Ay  3  8  *  10“5[m]  which  is  in  the  required  interval. 

For  both  the  laminar  and  turbulent  cases,  cells  in  the  z  direction  were 
adjusted  so  that  the  points  where  possible  physical  phenomena  such  as  shock 
waves  and  expansion  fans  are  expected,  very  fine  cells  were  used.  In  the 
other  parts  of  the  domain  larger  cells  were  used. 
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9 .  Heat  Transfer  Analysis 


Skin  friction  and  heat  transfer  quantities  were  calculated  in  both 
laminar  and  turbulent  cases  and  they  are  shown  in  Figures  (6  *  11). 

9 . 1  Laminar  Calculation 

In  laminar  flow  fluxes  can  be  derived  directly  from  the  gradients  near 
the  wall.  The  first  cell  is  close  "enough"  to  the  wall  and  gradients  of 
velocity  and  enthalpy  do  not  change  much  in  this  region  near  the  wall.  The 
shear  stress  and  heat  flux  in  the  laminar  case  will  be: 


W1 


Tw  "  U  AY]72 

(7.1) 

u  ^1  -  hw 
qw  "  77  AYj/2 

(7.2) 

The  skin  friction  coefficient  and  Stanton  number  will  be: 

r  =  2*tW 

Cf  p  U2 

(9.1) 

CO  oo 

(9.2) 

where  hr  is  the  recovery  enthalpy 

r(jzl)  2 
hr  14-  2 

ho  1  +  (JCI)  7 

(9.3) 

2  oo 

r-  is  the  recovery  factor 

r  -  /Pr  (laminar  flow) 

(9.4) 
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The  coefficient  of  heat  transfer  in  convection  was  calculated  using 


S  st 


(9.5) 


9 • 2  Turbulent  Calculations 

In  turbulent  flow  the  gradients  of  velocity  and  enthalpy  near  the  wall 
are  very  steep  and  change  rapidly  with  distance  from  the  wall. 

Direct  calculation  of  flux  gradients  is  not  accurate  in  this  case.  The 
log  law  approach  is  used  to  calculate  skin  friction.  In  the  calculations 
using  PHOENICS  flow  field,  the  following  relation  has  been  used. 


Cf 


2  k*, 


(9.6) 


w2P<»3 .33 

CO 

To  obtain  equation  9.6,  the  turbulent  kinetic  energy  equation  has  been  used  as 
a  starting  point.  [Ref.  5], 

Dk  (Vt  3k  f  7v 

P  Dt  3y  (Sk  3y  (9*7) 

+  k  [  (|isi>  .  cD  1 

k  v8y  u 

The  source  terra  of  the  turbulent  kinetic  energy  equation  should  be  zero  near 
the  wall  which  means 


1,1  '  3“  )?,  - 


k  ^  3y  'w 


(9.8) 


therefore  the  shear  stress  on  the  wall  can  be  defined  as: 

Tw  =  CD  172  Pw  kw  (9.9) 

where  is  the  turbulent  kinetic  energy  on  the  wall,  pw  is  the  density  on  the 
wall  and  Cq  *  0.09  [Ref.  5],  substituting  the  values  above  into  the  Blasius 
skin  friction  relation  the  Cf  equation  becomes: 
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2  Tw 


(9.10) 


Cf 


Poo 


00 


Pw  2  kw 


00 


The  heat  transfer  quantities  are  evaluated  from  the  Chil ton-Colburn  form  of 
Reynolds  analogy. 

St  -  (Cf / 2 )  *  Pr-2/3  (9.11) 

qw  -  st*pB*Ua*(hr-hw)  (9.12) 

where  equation  (9.3)  is  used  to  evaluate  hr  with  the  recovery  factor  given  as: 
1  /  3 

r  =■  Pr  (turbulent  flow)  (9.13) 

The  convective  heat  transfer  coefficient  is  calculated  by  using  equation  (9.5) 
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10.  Code  and  Computer 


PHOENICS  81,  Body  Fitted  Coordinate  (BFC)  version  has  been  used  in  the 
computations  (see  Ref.  3).  PHOENICS  has  been  installed  on  NPS  IBM  3033  MVS 
1.3  computer.  400  sweeps  per  computer  run  provided  a  reasonable  convergence 
in  all  runs  except  the  turbulent  blunt  case  continuity  error  of  less  than 
4.10”4  has  been  achieved  in  the  three  runs. 

The  continuity  error  is  the  total  summation  of  the  absolute  mass 
imbalance  in  all  cells  divided  by  the  inlet  mass  flux.  CPU  time  consumption 
varies  from  case  to  case  as  follows: 


Laminar  Wedge 

630 

CPU 

Seconds 

Turbulent  Wedge 

630 

CPU 

Seconds 

Laminar  Blunt 

630 

CPU 

Seconds 

Turbulent  Blunt 

1542 

CPU 

Seconds  for  1000  sweeps 

11.  Results  and  Discussion 


The  results  of  the  calculations  are  available  on  appendix  c.  The  tabular 
results  include  the  values  of  pressure,  velocities,  enthalpy,  temperature  mach 
number,  density,  turbulent  kinetic  energy  and  rate  of  turbulent  dissipation. 
The  values  are  given  in  18  x  29  cells  points. 

Skin  friction  and  heat  transfer  results  are  shown  in  Figures  (5-11). 
Laminar  and  turbulent  skin  friction  and  Stanton  number  in  wedge  flow  show 
improvement  compared  to  the  results  reported  by  Yukselen  (Ref.  10).  The  lines 
are  smoother  and  the  oscillations  at  the  end  were  eliminated.  Basically  the 
magnitudes  are  similar  to  those  in  Ref.  10. 

Laminar  blunt  values  are  similar  except  near  the  beginning.  The 
beginning,  as  expected  in  blunt  zone,  creates  higher  rates  of  heat  transfer. 
Even  though  the  blunt  geometry  used  is  a  multi-wedge  shape  it  should  predict 
the  correct  values  except  for  the  stagnation  point  itself. 

Turbulent  blunt  skin  friction  has  different  behavior.  It  has  a  very 
large  value  at  the  first  point  and  then  undershoots  to  values  that  are  smaller 
than  for  wedge.  It  should  also  be  kept  in  mind  that  the  convergence  of  this 
case  wasn't  very  successful. 
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12. 


rrmri nsions  and  Recommendations 


1.  PHOENICS  was  found  to  be  a  friendly  code  for  simulating  complicated 

mixed  heat  transfer  fluid  dynamics  problems. 

2.  Derivation  of  heat  transfer  properties  to  .a  vane  solid  wall  in 
laminar  and  turbulent  flow  has  been  installed  in  the  code.  It  can  be  used  for 
predictions  of  heat  transfer  rate  in  both  cold  and  hot  gas  flow. 

3.  Two  features  have  been  added  to  the  code  in  NPS:  The  restart  option 
and  the  use  of  initial  field,  make  it  possible  to  simulate  time  dependent 
processes  and  solve  the  temperature  variation  in  the  vane  itself. 
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Figure  1:  NWC  Jet  vane  configuration. 
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Figure  2-  NWC  Jet  Vane  Approximation 
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Figure  3:  Wedge  vane  domain  and  grid 
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Figure  4:  Blunt  vane  domain  and  grid. 
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Figure  5:  R  No.  along  the  Vane 
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Figure  6:  in  Laminar  flow. 


32 


A  C, 


io.io'34 


9 . 10'3+ 


8. 103+ 


7. 10"3+ 


6.10 


-3 


5.10 


-3 


-3 


4.-10 


3.10 


2.10 


-3 


1.10 


-3 


1 

2 


I 

3 


I 

4 


I  III  mill  III  I  I  l  I  .  '  r  7 

3  6  789  161718  19  20  21  22  23  27  28  29 

►  Z 


I 


I  I  1 1  I  I  I 


I 


- 1 - 1  II - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1  I  T 

0  10  20  30  40  50  60  70  80  90  100  110  120  130  140  150 

Figure  7-  in  Turbulent  flow 


33 


A  s 


t 


5 . 1 0”3+ 


M  =  3.2 

GO 

A  -  Wedge  Vane 
□  ~  Blunt  Vane 


I  I  I  I  I  II  I  IIUI  III  I  I  I  I  I  I  I  I  I  1  T  7 

1  2  3  4  5  6789  161718  19  20  21  22  23  2627  28  29  iL 

i  i - 1 - 1  I  l 1  I  "  l  I  i  i  »  i  i  ►  z 

0  10  20  30  40  50  60  70  80  90  100  110  120  130  140  (mm) 

Figure  8:  in  Laminar  flow 


34 


rtt^RODucCu  AT  GuVlRnmENT  EXPENSE 


8.IO'3- 
7.IO'3- 
6.10"3- 
5.1 0  3- 
4.10'3- 
3.1 0~3 
2.10 3 
1  .10‘3 


0 


0 


1  i  i  iin  nan  ill  i  i  i  1 1 1  i  i  i 

2  3  4  5  6  789  161718  19  2  0  21  22  23  2627  28 

T — ■T - T - 1 - 1  1  1  1  1  1  1  1  1  ^(mml 

10  20  30  40  50  60  TO  80  90  100  110  120  130  140 


Figure  9:  Turbulent  stonton  number. 


35 


t  h  l1"2*] 


8000- 


7000- 


6000- 


5000- 


4000- 


3000- 


2000 


1000 


0A 


r 

0 


I  I  I  I  II I  IIRRII  III  I  I 

Z  3  4  5  6  789  161718  19  20  21 


l  l  1 1  I  I  I  I 
22  23  27  28  29 


IZ 


i 

50 


T 


100 


>  z 

(mm) 


Figure  10  :  Coefficient  of  heat  convection  in  laminar  flow. 
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Appendix  A 


Satellite  Listing 

Two  subroutines  SATELLITE  and  GROUND  had  to  be  changed  and  improved.  The 
full  list  is  enclosed  in  Appendix  A  and  B. 

VAN4SAT  and  VANTSAT  are  the  laminar  and  turbulent  SATELLITE  subroutines, 
the  first  has  the  blunt  geometry  and  the  second  has  the  wedge  (it  can  be 
changed  easily  from  wedge  to  blunt  and  vice  versa)  VAN4GRD  and  VANTGRD  are 
exactly  the  same.  They  are  the  GROUND  subroutines,  VANTGRD  is  given  in 
Appendix  B. 
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file*  vantsat  fortran  Al 


C$ DIRECTIVExxSATL IT  AMI  LEITNER 

C  LAMINAR  SOLUTION  FOR  NWC5  NY=18  NZ=29  YN=GTH 

C  LECSAT  CONVERTED  TO  DIAMSAT 

C  xFILE  NAME*  MODBFCST . FTN 

C  ^ABSTRACT  *  SATELLITE  MODEL  MAIN  PROGRAM.  THIS  VERSION  IS 
C  FOR  USE  WITH  THE  BODY-FITTED  COORDINATE  SCHEME  (SUMMER  1984 
C  VERSION)  PROVIDED  AS  AN  ATTACHMENT  TO  SPRING  1983  PHOENICS. 

C  X DOCUMENTATION i  PHOENICS  INSTRUCTION  MANUAL  (SPRING  1983) 

C  NITH  BODY-FITTED  COORDINATES  INSTRUCTION  SUPPLEMENT 

C  (SUMMER  1984). 

C  XAUXILIARY  SUBROUTINES  (TAPES,  ETC.)  ARE  IN  SATELLITE  LIBRARY 

C  SERVICEU,  NHICH  MUST  BE  INCLUDED  IN  LINK  EDIT  TO  RUN. 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  1  STARTS: 


C - 

CHAPTER  1  COMMON  BLOCKS  AND  USER'S  DATA. 
C - - - 


INCLUDE  (CMNGUS) 

INCLUDE  (CMNGRF) 

INCLUDE  (GUSSEQ ) 

COMMON/CPI/ IPNRIT, I DUM( 243) 

DIMENSION  GDTAPE( 3) , DFAULT( 4) 

DIMENSION  ARRAY1 ( 309 ) , ARRAY2( 194) , ARRAY3( 421 ) 

LOGICAL  ARRAY1 , LSPDA, HRT, RD, NAMLST 
INTEGER  ARRAY2, XP LANE, YP LANE, ZP LANE 

INTEGER  P1,PP,U1,U2,V1,V2,W1,W2,R1,R2,RS,EP,H1,H2,H3,C1,C2, 
&C3, C4 

REAL  NORTH, LON 
LOGICAL  BFC 

EQUIVALENCE  ( ARRAY1 ( 1 ) , CARTES ) , ( ARRAY2( 1 ) , NX ) 

EQUIVALENCE  ( ARRAY3( 1 ) , SPARE1 ( 1 ) ) , ( Ml , R1 ) , ( M2 , R2 ) 
EQUIVALENCE  ( LSTRUN, INTGR( 12) ),( NAMLST, L0GIC(88) ) 
EQUIVALENCE  ( LOGIC( 20 ) , BFC) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  1  ENDS. 
C$DIRECTIVE**CMNBF1$$ 

C  THIS  FILE  CONTAINS  SATELLITE  COMMON  BLOCKS  FOR  BFC'S 
C  FI  MUST  BE  DIMENSIONED  TO  GREATER  THAN  OR  EQUAL  TO 

C  (NX+NY+17xNZ+24*NX*NY+6X(NX+1)x(NY+1)+6xND) *  THE  VALuE 

C  OF  THE  DIMENSION  MUST  BE  SET  AS  NBFC  IN  GROUP  6  OF  SATLIT. 
COMMON/ FOB/ FI ( 5000 ) 

COMMON/CIB/ND/CIC/KOORD 

COMMON/CID/KDBGG,KDBGMF,KDBGCD,KDBIND,KDBMFX,KDBCDT,KDBPCS, 
&  KDBGUV, KDBGPV 

COMMON/CI E/KDBGS, KDBINS 
COMMON/ Cl F/IGEN/CIG/NCART 

C  THE  FOLLONING  ARRAYS  MUST  BE  EXACTLY  DIMENSIONED  FOR  NXP1, 

C  NYP1  AND  NZP1 ,  BUT  MAY  BE  OVER  DIMENSIONED  FOR  ND. 

C  THE  BFRAC  ARRAYS  MUST  BE  DIMENSIONED  TO  ALLOW  FOR  SETTINGS 
C  IN  SATLIT,  THEY  MAY  BE  OVER  DIMENSIONED. 

COMMON/CRA/XW( 19,30,1 )/CRB/XE( 19 , 30,1) 

&  /CRC/YS (2,30, 1)/CRD/YN(2,30, 1) 

&  /CRE/ZL ( 2, 19 , 1 )/CRF/ZH( 2, 19 , 1 ) 

&  /CRG/RC0N/CRH/DARCY/CRI/BXFRAC(99)/CRJ/BYFRAC(99) 

&  /CRK/BZFRAC(99) 

COMMON/CL A/STORSA( 6 ),STORWD( 6 ) , STORP, STORPE, STORPN, 

&  STORPH, STOR1 , ST0R2, ST0R3, STOUNV, PRTBFC, STOCRN 

COMMON/CLC/BFPLOT 

LOGICAL  STORP, STORPE, STORPN, STORPH, STOR1 , ST0R2, ST0R3, 

&  STORSA, STORWD, STOUNV, PRTBFC, BFPLOT , STOCRN 


C  END 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  USER  SECTION  1  STARTS* 

C  GRAFFIC  ARRAYS  DIMENSIONED  AS  NEEDED... 

COMMON/GRAF1/PHI1 ( 1 )  /GRAF2/PHI2( 1 ) 

C  .  POROSITY  &  SPECIAL  DATA  ARRAYS  DIMENSIONED  AS  NEEDED... 
DIMENSION  PE(1,1,1),PN(1, 1,1), PH (1,1,1), PC (1,1,1) 

DIMENSION  LSPDA(l) , ISPDA( 1 ) , RSPDA( 1 ) 

C  USER  PLACES  HIS  VARIABLES,  ARRAYS,  EQUIVALENCES  ETC.  HERE. 
EQUI VAL  ENCE( RAIR,RE(21)), (GAMA, RE( 22) ) , (GSWP, RE( 23 ) ) 

1, (GPR,RE(24)),(TW,RE(25)), (GEMU1 , RE(26 ) ) , ( JEMU1 , INTGR( 1 ) ) 

C  USER  PLACES  HIS  DATA  STATEMENTS  HERE. 

DATA  NLSP,NISP,NRSP/ 1,1,1/ 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  USER  SECTION  1  ENDS. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  2  STARTS* 


VAN0001 0 
VAN00020 
VAN00030 
VAN00040 
VAN00050 
VANOO  060 
VAN00070 
VAN0008  0 
VAN00090 
VAN00100 
VANOO 1 10 
VAN00120 
VAN00130 
VAN00140 
VAN00150 
VAN00160 
VAN0017  0 
VAN00180 
VAN00190 
VAN00200 
VAN0021 0 
VAN00220 
VANO  0230 
VAN00240 
VAN00250 
VANO  0260 
VAN0027  0 
VAN00280 
VANO  029  0 
VAN00300 
VAN00310 
VAN00320 
VAN00330 
VAN00340 
VANOO  350 
VANOO  36  0 
VAN0037  0 
VAN00380 
VANO  039  0 
VAN00400 
VAN0041 0 
VAN00420 
VAN00430 
VANO  0440 
VAN00450 
VAN00460 
VAN0047  0 
VAN00480 
VAN00490 
VANO  050  0 
VAN00510 
VAN00520 
VAN00530 
VAN00540 
VANO  0550 
VAN00560 
VAN00570 
VAN00580 
VANO  059  0 
VAN00600 
VAN00610 
VAN00620 
VAN00630 
VAN00640 
VAN00650 
VAN00660 
VAN0067  0 
VAN00680 
VAN00690 
VAN007  00 
VAN007 1 0 
VAN00720 
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FILE*  VANTSAT  FORTRAN  A1 


C - 

CHAPTER  2  SET  CONSTANTS,  AND  ARRANGE  FILE  MANIPULATIONS. 

C  PLEASE  DO  NOT  ALTER,  OR  RE-SET,  ANY  OF  THE  REMAINING 
C  STATEMENTS  OF  THIS  CHAPTER. 

DATA  CELL, EAST, WEST, NORTH, SOUTH, HIGH, LOW, VOLUME/ 

8  0.,1.,2.,3.,4.,5.,6.,7.  / 

DATA  P1,PP,U1,U2,V1,V2,W1,W2,RI,R2,RS,KE,EP,HI,H2,H3,C1,C2, 
8C3 , C4/1 ,2,3,4 ,5,6,7 ,8, 9, 10, 11, 12, 13, 14, 15, 16, 17 , 18, 19 ,20/ 

DATA  FIXFLU, FIXVAL , ONLYMS, WALL/1. E-l 0,1 .El 0,0 .0,-10 . 0/ 

DATA  IPLANE, XPLANE, YPLANE, ZPLANE/0 , 1,2,3/ 

DATA  WRT, RD, DFAULT/ .TRUE .,. FALSE . , 4HDEFA, 4HULT . , 4HDTA/, 1HG/ 
DATA  GDTAPE/4HGUSI , 4HE1 . D, 2HTA/ 

DATA  NLDATA, NIDATA, NRDATA/309 , 194, 421/ 

DATA  NLCREG, NTCVRG/60 ,350/ 

DATA  TITPP,TITC1, TITC2, TITC3/3HRH0, 4HMACH, 4HTEMP, 4HCFST/ 

CALL  TAPES(10,GDTAPE,3,1,4*NRDATA) 

c - READ  DEFAULT  FILE  IF  BLOCKDATA  ABSENT 

IF(INTGR1(29).NE.10)  GO  TO  2 

CALL  WRIT40C 40HDATA  ESTABLISHED  IN  BLOCK  DATA.  ) 

GO  TO  3 

2  CALL  DEFLT 

CD  2  CALL  TAPESC 1 , DFAULT , 4 , 2, 4*NRDATA) 

CD  CALL  DATAIOC  RD, 1 ) 

CALL  WRIT40C40HDATA  TAKEN  FROM  DEFAULT. DTA  ON  GROUP  A/C) 

3  CALL  WRIT40C40HFILE  MODSTL.FTN  IS  THE  SATLIT  USED.  ) 

L0GICC89 X . TRUE . 

C - 

CHAPTER  3  DEFINE  DATA  FOR  NRUN  RUNS. 

C - 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  2  ENDS. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  USER  SECTION  2  STARTS* 

C—  GROUP  41MULTI-RUNS  *  RUN(  1-30  X  .T  .  ,  29* .  F .  > 

C 

RUNC1X. FALSE. 

C  NOTE*  ALL  RUNS  ARE  DEACTIVATED  AT  THIS  POINT  -  USER  SHOULD 
C  ====  SWITCH  ON  ONE  ONLY  OF  RUNS  1-4  IN  NEXT  STATEMENT. 

RUN( 1 X . TRUE . 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  USER  SECTION  2  ENDS. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  3  STARTS* 
DO  10  IRUN=1,30 

I FC .NOT.RUN(IRUN))  GO  TO  10 
NRUN=NRUN+1 
LSTRUN=IRUN 
10  CONTINUE 

DO  999  IRUN=1 , LSTRUN 

IF( .NOT.RUN(IRUN))  GO  TO  999 
INTGR(ll)  =  IRUN 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  3  ENDS. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  USER  SECTION  3  STARTS* 

C—  ALL  INTEGER  VARIABLES  ARE  DEFAULTED  TO  0,  AND  REAL  VARIABLES 
C  TO  0.0,  UNLESS  OTHERWISE  INDICATED. 

C  E.G.  BY  VARIABLE<10>,  OR  <10. 0>  AS  APPROPRIATE. 

C  THE  DEFAULT  SETTINGS  OF  ALL  LOGICAL  VARIABLES  ARE  ALWAYS 
C  INDICATED,  E.G.  VARIABLE< . T . >,  OR  VARIABL E< . F . > . 

C 

C -  RUN1 

C -  -  - 

C—  GROUP  1.  FLOW  TYPE  * 

C  PARAB< . F. >, CARTES< . T . > , ONEPHS< . T . > 

C - -  -  - 

C—  GROUP  2.  TRANSIENCE  * 

C  STEADY<  . T  .  >,  ATIME,  LSTEP<1>,  FSTEP<1> 

C  *  TLAST<1 . E10> , TFRACC 1-30  X30X1  .  > 

C  SERVICE  SUBROUTINE  FOR  1  NT  *  POWER-LAW  TIME  STEPS* 

C  CALL  GRDPWRCO, NT, TLAST, POWER) 

C - -  -  - - 

C—  GROUP  3.  X-DIRECTION  * 

C  NX<1>,XULAST<1 . 0>,XFRAC( 1-30) 

C  SERVICE  SUBROUTINE  FOR  POWER-LAW  GRID* 

C  CALL  GRDPWRC1, NX, XULAST, POWER) 

C -  -  - 


VAN00730 
VAN00740 
VAN00750 
VAN00760 
VAN0077  0 
VAN00780 
VAN00790 
VAN00800 
VAN00810 
VAN00820 
VAN00830 
VAN00840 
VAN00850 
VAN00860 
VAN00870 
VAN00880 
VAN00890 
VAN00900 
VAN00910 
VAN00920 
VAN00930 
VAN00940 
VAN00950 
VAN0096  0 
VAN0097  0 
VAN00980 
VAN00990 
VAN01000 
VAN01010 
VAN01020 
VANO 1030 
VAN01040 
VAN01050 
VAN01060 
VAN0107  0 
VANO 1080 
VAN01090 
VAN011 00 
VAN0111 0 
VANO 1120 
VAN01130 
VANO 1 140 
VANO 1 150 
VANO 1 16  0 
VANO 1 17  0 
VAN01180 
VANO 1 190 
VANO 1200 
VAN01210 
VAN01220 
VAN01230 
VAN01240 
VAN01250 
VANO 1260 
VANO 1270 
VAN01280 
VAN01290 
VANO 1 300 
VAN01310 
VAN01320 
VAN01330 
VAN01340 
VAN01350 
VAN0136  0 
VAN0137  0 
VAN01380 
VANO 1390 
VANO 1400 
VAN01410 
VAN01420 
VANO 1430 
VAN01440 
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FILE*  VANTSAT 


FORTRAN  A1 


-  GROUP  4.  Y-DIRECTION  i 

NY<1>,YVLAST<1 .0>,YFRAC(1-30),RINNER,SNALFA 
SERVICE  SUBROUTINE  FOR  POWER-LAW  GRID* 

CALL  GRDPWRC2, NY, YVLAST, POWER) 

NY=18 


-  GROUP  5.  Z-DIRECTIQN  : 

NZ<1>,ZWLAST<1 . 0>, ZFRAC(l-30 ) 

SERVICE  SUBROUTINE  FOR  POWER-LAW  GRID* 
CALL  GRDPWR(3,NZ,ZWLAST, POWER) 

NZ=29 


-  GROUP  6.  MOVING  GRID  OR  DISTORTED  (BODY-FITTED*)  GRID  : 

—  MOVING  GRID  : 

MGRID, IZW1 , IZW2, AZW2, BZW2, CZW2, PINT#  ZW2M1T 


-  BODY-FITTED  GRID  — 

BFC< .T. >, IGEN<1>,  ND<1>,NBFC<5000>,KQQRD,RCON 
BXFRACC  l-NXXl .  0 , NXMlxO  .  0> 

BYFRACd-NYXl  .0,NYM1*0 . 0> 

BZFRACC1-NZX1 .0,NZM1*0 . 0> 

SERVICE  SUBROUTINE  FOR  SUB-DOMAIN  SPECIFICATION  (FOR  IGEN=1 
ONLY) : 

CALL  DOMAIN(ID,IXF,IXL,IYF,IYL,IZF,IZL) 

XE(1-NYP1,1-NZP1, 1-NDX(NYP1*NZPI*ND)x1 .  0>, 

XW( 1-NYP1 , 1-NZP1 , 1-ND) , 

YN(  1-NXP1  ,  1-NZP1 , 1-NDX(NXP1*NZP1XND)X1 .  0>, 
YS(1-NXP1,1-NZP1,1-ND), 

ZHd-NXPl,  1-NYP1,  l-NDX(NXPlxNYPlxND)*l  .0>, 

ZL(1-NXPI,1-NYP1,  1-ND),5TQRSA(  1-6  X6*.F.>,  STQRWD(  1-6X6  x.F.>, 
STORP<.F.>,STQRPE<.F.>,STQRPN<.F.>,STQRPH<.F.>,STOUNV<.F.>, 
PRTBFC< . F. >, DARCY, BFPLOT< . F. > 

CYCLIC  BOUNDARY  CONDITIONS  ARE  DEFAULTED  INACTIVE  ; 

TO  ACTIVATE  THEM  AT  SELECTED  IZ  SLABS  USE  SERVICE  SUBROUTINE : 
CALL  XCYIZ(IZ, .TRUE. ) 

SERVICE  SUBROUTINE  TO  DEACTIVATE  CURVATURE  TERMS  IN  U,  V 
AND  W  EQUATIONS  ASSOCIATED  WITH  CURVATURE  OF  IX,  IY,  IZ 
GRID  LINES  RESPECTIVELY: 

CALL  UCURVEdZ,  .FALSE.  ) 

CALL  VCURVEUZ,  .FALSE.  ) 

CALL  WCURVEdZ,  .FALSE.  ) 

NCART<1> 

^WARNINGS  I  I  |  I  I  I 


A)  WHEN  USING  BFCS  ST0VAR(H3),  STQVAR(C4),  ST0VAR(21)  ARE 
AVAILABLE  ONLY  FOR  STORING  NQN-QRTHQGONAL  VELOCITY 
COMPONENTS. 

B)  MULTI-RUNS  ARE  NOT  ALLOWED  WITH  BFC  OPTION. 

C)  MOVING  GRID, TWO-PHASE  AND  PARABOLIC  OPTIONS  ARE  NOT 
AVAILABLE  WITH  BFC  OPTION. 

D)  KE-EP  TURBULENCE  MODEL  SHOULD  BE  USED  WITH  BFC'S  ONLY 
WHEN  THE  MAIN  FLOW  IS  IN  THE  IZ  DIRECTION. 

E)  BUILT-IN  GRAVITY  TERMS  DO  NOT  TAKE  ACCOUNT  OF  BFC'S. 
XNOTES 


A)  THE  STANDARD  VELOCITY-FIELD  PRINTOUT  FOR  THE 
VELOCITY  RESOLUTES  IS  ACTIVATED  IN  THE  USUAL 
WAY.  AN  ADDITIONAL  OPTION  EXISTS  FOR  PRINTING  THE 
CARTESIAN  VELOCITY-COMPONENTS  WHICH  MAY  BE 
ACTIVATED  BY  SETTING  THE  FOLLOWING  LOGICALS: 

STQVAR(U2)=.T.  FOR  U-COMPONENT  (CARTESIAN) 
ST0VAR(V2)=.T.  FOR  V-COMPQNENT  (CARTESIAN) 

STQVAR( W2)= . T .  FOR  W-COMPONENT  (CARTESIAN) 

SIMILARLY  PRINTOUT  OF  NON-QRTHQGGNAL  VELOCITY 
COMPONENTS  MAY  BE  ACTIVATED  AS  FOLLOWS: 

STQVAR(C4)=.T.  FOR  U-COMPONENT  (NQN-ORTHOG) 
STQVAR(H3)=.T.  FOR  V-COMPONENT  (NON-ORTHOG) 

STQVAR( 21 ) = . T .  FOR  W-COMPONENT  (NON-ORTHOG) 

B)  BFC  (TO  ACTIVATE  THE  BFC  OPTION),  IGEN  (THE  CODE  FOR  METHOD 

OF  GRID  SPECIFICATION),  ND  (NUMBER  OF  SUB-DOMAINS)  AND 
NBFC  (THE  FI  ARRAY  DIMENSION),  MUST  BE  SET  BEFORE 
"STANDARD  BFC  SECTION  2".  ====== 


VAN01450 
VAN0146  0 
VAN0147  0 
VAN01480 
VAN01490 
VANO 1500 
VAN01510 
VANO 1520 
VAN01530 
VANO 1540 
VAN01550 
VAN0156  0 
VAN0157  0 
VAN01580 
VAN01590 
VAN016  00 
VANO 16 10 
VANO 1620 
VAN016  30 
VANO 1640 
VAN01650 
VANO 166  0 
VAN01670 
VAN01680 
VAN01690 
VAN017  00 
VANO 1710 
VANO 1720 
VAN017  30 
VANO 1740 
VANO 1750 
VAN01760 
VANO 177  0 
VAN01780 
VANO 1790 
VANO 1800 
VAN0181 0 
VANO 1820 
VAN01830 
VANO 18  40 
VANO 1850 
VANO 186  0 
VANO 187  0 
VAN01880 
VAN0189  0 
VANO 190  0 
VAN0191 0 
VANO 1 920 
VAN01930 
VAN01940 
VAN01950 
VANO 196  0 
VANO 197  0 
VANO 1980 
VANO 1 990 
VAN02000 
VAN02010 
VAN02020 
VAN02030 
VAN02040 
VAN02050 
VAN02060 
VAN02070 
VAN02080 
VAN02090 
VAN02100 
VAN02110 
VAN02120 
VAN021 30* 
VAN02140 
VAN02150 
VAN0216  0 
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FILE*  VANTSAT  FORTRAN  A1 


C  ALL  OTHER  BFC  DATA  MUST  BE  SET  AFTER  "STANDARD  BFC 

C  SECTION  2.  ===== 

C  C)  NXP1,  NYP1 ,  NZP1  STORE  NX+1,  NY+1,  NZ+1;  THESE  ARE 

C  AVAILABLE  TO  USER  AFTER  STANDARD  BFC  SECTION  2. 

C  D)  FOR  IGEN=1  USE  BXFRAC, BYFRAC  &  BZFRAC  IN  PLACE  OF 

C  XFRAC/YFRAC  &  ZFRAC . 

C - 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  BFC  SECTION  1  STARTS i 
C  DEFAULT  SETTINGS i 
NCART=10 
BFC= . TRUE . 

IGEN=1 

ND=1 

NBFC=5000 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  BFC  SECTION  1  ENDS.. 

C  XUSER  SETS  BFC,  IGEN,  ND  AND  NBFC  HERE. 

C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  BFC  SECTION  2  STARTS. 
CALL  SB4I C NXP1 , NX+1 , NYP1 , NY+1 , NZP1 , NZ+1 , 1,0) 

IFCBFC)  CALL  BFCDFTCNBFC,XE,XW,YN,YS,ZH,ZL,ND,NXP1,NYP1, 

8  NZPl.NZ) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  BFC  SECTION  2  ENDS. 

C  XUSER  SETS  ALL  OTHER  BFC  VARIABLES  HERE. 

C  XUSING  NONIFORM  GRID  1-8 
GTH=65.E-3 
GTL=150 . E-3 
GBETA=4. 

GBETA=GBETAX3. 1415927/180 
GTAB=TAN(GBETA) 

DELMAX=2 . E-3 
GNBL=5 . 

GPWR=2 . 

DO  64  IY=1 , 5 

64  BYFRAC(IY)=CFLOAT(IY)/GNBL)XXGPWRXDELMAX/GTH 
BYFRACC6 )=BYFRACC5)+3 . E-3/GTH 

DEL  =  C 1 .-BYFRACC6 ) )/( FLOAT  C  NY)-GNBL-1 ) 

DO  65  IY=7 ,NY 

65  BYFRAC( IY) =BYFRACC IY-1 )+DEL 

BZFRACC 1 ) =10 . E-3 
DO  66  IZ=2,5 

66  BZFRACC IZ) =10. E-3+BZFRACC IZ-1) 

BZFRAC(6)=BZFRACC5)+5.E-3 

DO  67  IZ=7 , 9 

67  BZFRAC(IZ)=BZFRAC(IZ-l)+2.E-3 
DO  68  IZ=10,10 

68  BZFRACCIZ)=BZFRACCIZ-l)+l.E-3 
DO  77  IZ=11 , 14 

77  BZFRACC I Z) = BZFRACC IZ-1 )+ .  5E-3 
DO  78  IZ=15, 15 

78  BZFRACC IZ) = BZFRAC ( IZ-1 )+l . E-3 
BZFRACC 16 ) =BZFRACC 15 )+2 . E-3 
BZFRACC 17 ) =BZFRACC 16 )+3 . E-3 
BZFRACC 18 )=BZFRACC 17 )+5. E-3 
DO  69  IZ=19, 22 

69  BZFRACCIZ)=BZFRACCIZ-l)+10.E-3 
BZFRACC  23 )= BZFRACC  22) +3 . E-3 
BZFRACC  24 ) =BZFRACC  23 )+2 . E-3 
BZFRACC 25) =BZFRACC24)+2. E-3 
BZFRACC 26 )=BZFRACC 25 )+3. E-3 
BZFRACC  27 ) =BZFRACC  26 )+5 . E-3 

DO  71  IZ=28 , NZ 

71  BZFRACC IZ)=BZFRACC IZ-1 )+10 . E-3 
DO  72  IZ=1,NZ 

72  BZFRACC IZ)=BZFRACCIZ)/GTL 
CALL  DOMAINC 1 , 1 , NX, 1 , NY, 1 , NZ) 

DO  61  IX=1 , NXP1 

DO  62  IY  =  1 , NYP1 
ZL  C IX , I Y, 1 ) =0 . 0 
62  ZHCIX, IY, 1)=GTL 
DO  63  IZ=1, NZP1 
YNCIX,IZ,1)=GTH 
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63  YS(IX,IZ,1)=0.0 

C  YS(IX,13,1)  SHOULD  COME  AFTER 
DO  662  IZ=16,25 

°662  YS(IX?Izfi)*(BZFRAC(IZ-l)-BZFRAC(3))»OTABXOTL 
DO  663  IZ=13,15 

GZ12=( BZFRAC( IZ-1 )-BZFRAC( 11 ) )*GTL 
663  YSCIX,IZ,1)=SQRTCYSCIX,16,1 )*GZ12#2 • "GZ12XX2) 


DO  664  IZ=26  , NZ 
664  YSCIX,IZ,1)=YSCIX,25, 1) 

61  CONTINUE 

STORSAC IFIX( LOW) )= .TRUE. 
STORSAC I FIXC HIGH) )= .TRUE. 
STORSAC IFIXC SOUTH) )=. TRUE. 
STORWDC IFIXC SOUTH) )=. TRUE. 
STORP= .TRUE . 

PRTBFC= . TRUE . 

CDAR  DARCY=1 . E10 


-  GROUP  7.  BLOCKAGE:  BL0CK< . F . >, IPLANE, IPWRIT 

XSET  CONSTANT  POROSITIES  OVER  SUB-DOMAINS  USING: 

CALL  CONPORCIR,TYPE,VALUE,IXF,IXL,IYF,IYL,IZF,IZL),  WHERE: 

IR  =  RUN  SECTION  NUMBER,  E.G.  1  FOR  RUN1  SECTION;  1  TYPE*  =  EAST, 
WEST,  NORTH,  SOUTH,  HIGH,  LOW  8  CELL.  1  VALUE  * =WANTED  POROSITY 
OVER  REGION  IXF, . . .IZL. 

^DIMENSION  ARRAYS  PEC  NX, NY, NZ) ,  PNC  NX, NY, NZ) ,  PHCNX, NY, NZ) ,  & 
PCC NX, NY, NZ)  ABOVE. 

XFOR  FULLY-BLOCKED  CELLS  CIE.  1  VALUE' =  0.0)  USER  NEED  SET  ONLY 
THE  'CELL '  POROSITY  CTO  ZERO),  AS  CELL-FACE  AREAS  ARE  THEN 
AUTOMATICALLY  ZEROED. 

XFOR  SATELLITE  PRINTOUT  OF  ALL  POROSITIES  IN  DOMAIN,  'IPLANE'= 
XPLANE  YPLANE  OR  ZPLANE,  FOR  DESIRED  CROSS-SECTION  DIRECTION. 
XFOR  EACH  'TYPE'  A  MAXIMUM  OF  10  CALLS  TO  CONPOR  IS  ALLOWED, 
BUT  IF  REQUIREMENTS  EXCEED  THIS  PROVISION  SET  BLOCK=.T.  8 
IPWRIT  =  -l ,  AND  SET  POROSITY  ARRAYS  EXPLICITLY  HERE  AS  WANTED. 
IN  THIS  CASE,  THE  USER  MUST  SET  A  L  L  ELEMENTS  OF 
ARRAYS  PE,  PN,  PH,  PC  C MANY  MAY  BE  0.0  OR  1.0).  HE  MAY  USE: 
CALL  CRCPARRAY, VALUE, IXF, IXL , IYF, IYL , IZF, IZL , NX, NY, NZ) 

ANY  NUMBER  OF  TIMES,  TO  SET  'PARRAY'  C=  PE,  ETC.)  TO 
'VALUE'  OVER  RANGE  IXF  TO  IXL,  IYF  TO  IYL,  IZF  TO  IZL. 

XCONPOR  MUST  NOT  BE  USED  IN  CONJUNCTION  WITH  EXPLICIT 
SETTINGS  OF  THE  ARRAYS  CINCLUDING  SETTINGS  VIA  CR). 


GROUP  8. DEPENDENT  VARIABLES  TO  BE  SOLVED  FOR  OR  STORED  : 
S0LVARC1-25X25X. F.>,STOVARC 1-25) <25X.F.>,C0NC1C 1-4 )<4X.T.> 
USE  FOLLOWING  NAMED  INTEGERS  FOR  ARRAY  ELEMENTS  1-20: 
P1,PP,U1,U2,V1,V2,W1,W2,M1,M2,RS,KE, EP, HI , H2, H3, Cl , C2, C3, C4 . 
S0LVARCP1)=.TRUE. 

SOLVARCPP)= .TRUE. 

S0LVARCV1)=.TRUE. 

SOL VARCW1)=. TRUE. 

S0LVARCH1)=.TRUE. 

SOL VARCKE)=. TRUE. 

SOLVARCEP)=.TRUE. 

STOVARC V2 ) = . TRUE . 

ST0VARCW2)=.TRUE. 

ST0VARCC1)=. TRUE . 

ST0VARCC2)=.TRUE. 

ST0VARCC3)=.TRUE. 


-  GROUP  9.  VARIABLE  LABELS  : 

TITLEC1-25)<2HP1,2HPP,2HU1,2HU2,2HV1,2HV2,2HW1, 2HW2 , 2HR1 , 
2HR2,2HRS,2HKE,2HEP,2HH1,2HH2,2HH3,2HC1, 2HC2, 

■  2HC3 , 2HC4 , 2HRX, 2HRY, 2HRZ,  2X4HXXXX> 

TITLECC1 ) =TITC1 
TITLECC2)=TITC2 
TITLECC3)=TITC3 
TITLECPP)=TITPP 


-  GROUP  10  PROPERTIES: 

IRH01<1>, IRH02<1>,RH01<1 . 0>,RH02<1 . 0>, 
ARH01<1 .0>,BRH01<1 .0>,CRH01<1 .0> 
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IEMU1<1>,EMUX1.0>,EMULAM<1.E-10> 

IHSAT, H1SAT , H2SAT , PSATEX<1 . 0> 
SIGMAC1-25X1.0,2.0,1.,1.E10,1.»1.E10,1.,1.E10 
4X1.0,1.314,1.0, 1.E10, 10X1. 0> 


IRH01=-1 
PT0T=55.E5 
T0T=555 . 55 
RAIR=287. 

GAMA=1 . 35 

CP=RAIR/( 1-1/GAMA ) 

TW=323 . 

HWALL=TWXCP 
HT0T=CPXT0T 
RHTOT =PT0T/T0T/RAIR 
L0GIC(87)=.TRUE. 

ARH01 =RHT0T/PT0TXX( 1/GAMA ) 
BRH01=1 ./GAMA 
TURBULENT  OR  LAMINAR 
IEMU1=2 


IEMU1=-1 
JEMU1=IEMU1 
EMU1=1 . E-5 
EMULAM=EMU1 
GEMU1=EMU1 
GPR= . 7 

SIGMA(24)=GPR 
SIGMAC 14  )  =  .  9 


—  GROUP  11  INTER-PHASE  TRANSFER  PROCESSES  i 

ICFIP,CFIPS, IMDOT ,CMD0T ,CA1I<1 . E6>,CA2I<1 . E6> 


-  GROUP  12  SPECIAL  SOURCES  > 

I SPCSO( 1 -25 ) , AGRAVX, AGRAVY, AGRAVZ, ABUOY, HREF 


-  GROUP  13  INITIAL  FIELDS  i 

FI INIT( 1-25) <25X1. E-10> 

MACH  NO.  OF  FREE  STREAM 
GMACH=3 . 2 

A=1+(GAMA-1 J/2XGMACHXX2 
TE=TOT/A 

RHE=RHTOT/Axx( 1/(GAMA-1 ) ) 

PSTAT=PTOT/AXX(GAMA/(GAMA-l)) 

RH01=ARH01XPSTATXXBRH01 

SONIC=SQRT(GAMAXRAIRXTE) 

WIN=SONICXGMACH 

RKEIN=0 .01XWINXX2 

EPIN=0 . 16XRKEINXX1 . 5/GTH/2 . 

FIINIT(W1)=WIN 

FIINIT(P1)=PSTAT 

FIINIT(Hl)=HTOT 

FIINIT(KE)=RKEIN 

FIINIT(EP)=EPIN 


-  GROUP  14  BOUNDARY/INTERNAL  CONDITIONS  : 

ILOOP1, IL00PN»XCYCLE< . F. >, P BAR, REGION! 1-10  XlOX . T . > 

XN.B.  ALL  10  REGIONS  ARE  DEFAULTED  .TRUE..  THE  USER  SHOULD 
SET  REGION(I)=. FALSE.  FOR  UNUSED  REGIONS  *1*. 

DO  14  1=1,10 
14  REGION(I)=. FALSE. 


—  GROUP  15  TO  24;  REGIONS  1  TO  10 

—  ONLY  THOSE  REGIONS  ARE  ACTIVE  WHICH  ARE  SPECIFIED  BY  THE 
.  USER,  PREFERABLY  BY  WAY  OFi- 

CALL  PLACE! I REGN, TYPE, IXF,IXL,IYF,IYL,IZF,IZL)  8 
CALL  COVAL( IREGN, VARBLE,COEFF, VALUE) 

CALL  PLACE(l,LOW,l,NX,l,NY,l,l) 

C  CALL  C0VAL(1,M1,FIXFLU,WINXRHE) 

CDAR  CALL  COVAL( 1 ,M1 , 1 . E-20, 1 . E+20XWINXRHE) 

GCM=2XGAMA/WI N/ ( GAMA-1 ) 

GVM=PTOTXRHE/ RHTOT 
CALL  C0VAL(1,M1,GCM,GVM) 

CALL  COVAL ( 1 , HI , ONLYMS, WIN) 
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CALL  C0VAL(1,H1,0NLYMS,HT0T) 

CALL  C0VAL(1,KE,0NLYMS,RKEIN) 

CALL  C0VAL(1,EP,0NLYMS,EPIN) 

CALL  PLACE(2,HIGH, 1,NX, 1,NY,NZ,NZ) 

C  CALL  COVAL (2,M1>FIXVAL,PSTAT*0. ) 

CALL  COVAL(2,M1,1000XWIN*RHE/PSTAT,PSTAT) 
CALL  C0VAL(2,H1,0NLYMS,HT0T) 

C  WALL  ALONG  THE  VANE  IZ(11,NZ) 

GCM=EMUl/( . 5XBYFRACC 1 XGTH) 
DYl=BYFRAC(l)xGTH 
GOEFF=EMUl/(0.5xDYl) 

GOEFH=EMU1/(0 .5XDY1XSIGMAC24) ) 

CALL  PLACE(3,S0UTH,1,NX,1,1,12,NZ) 

CALL  COVAL (3,W1,GOEFF,0.) 

CALL  C0VAL(3,H1,G0EFH,HWALL) 

CALL  COVAL (3,W1,WALL,0. ) 

CALL  COVAL ( 3, HI , WALL ,  HWALL  X 
CALL  COVAL(3,KE,WALL,0. ) 

CALL  COVAL(3,EP,WALL,0. ) 


—  GROUP  25  GROUND  STATION  i 
GR0STA< . F . > , NAMLST< .  F .  > 

XNAMLST  ACTIVATES  NAMELIST  IN  GROUND. 
GR0STA= .TRUE. 


-  GROUP  26  SOLUTION  TYPE  AND  RELATED  PARAMETERS  * 

WHOL  EP< . F . >, SUBPST< . F . > , D0NACC< . F . > 

WH0LEP= . TRUE. 


-  GROUP  27  SWEEP  AND  ITERATION  NUMBERS  i 

FSWEEP<1>,LSWEEP<1>,LITHYD<1>,LITC<1>,LITKE<1>,LITH<1>, 
LIT  ERCl -25X9X1, -1,15X1> 

IVELF<1>,NVEL<1>,IVELL<10000>, 

IKEF<1>,NKE<1>, IKEL<1 0000> , 

IENTF<1>,NENT<1>, I ENTL<1 0000>, 
ICNCF<1>,NCNC<1>,ICNCL<10000>, 
IRH01F<1>,NRH01<1>,IRH01L<10000>, 

IRH02F<1>,NRH02<1>, IRH02L<10000> 

LSWEEP=1201 
GSWP=LSWEEP 
FSWEEP=801 
L ITERC  PP) =20 
L ITER( VI ) =5 
LITERCWl )=5 
LITHYD=2 


-  GROUP  28  TERMINATION  CRITERIA  * 

ENDITC  1-25X9x1 . E-10, 0 . 5, 15X1 . E-10> 
ENDIT ( 1 ) =1 . E-5 


-  GROUP  29  RELAXATION  : 

RLXP<1 . >,  RLXPXY<1 . > p RLXPZ<1 . >, RLXRH0<1 . >, RLXMDT<1 . >, 
DTFALSC3-25X23X1  .E10> 

DTFALS(Wl)=l.E-5 
DTFALSC VI )=1 . E-5 
DTFALS(KE)=1 .E-5 
DTFALSC  EP) =1 . E-6 
RLXP= . 3 


-  GROUP  30  LIMITS  i 

VELMAX<1 .E10>,VELMIN<-1 . El 0>, RH0MAX<1 . El 0>, RH0MIN<1 . E-10>, 
TKEMAX<1 . E10>, TKEMIN<1 . E-10>, EMUMAX<1 . E10>, EMUMIN<1 . E-10>, 
EPSMAX<1 . E10>, EPSMIN<1 . E-l 0>, AMDTMX<1 . E10>, AMDTMN<-1 . E10> 
EPSMAX  =  1 . El  3 


-  GROUP  31  SLOWING  DEVICES  :  SL0RH0<1 . > , SL0EMU<1 . > 

SL0RH0= . 2 


-  GROUP  32  PRINT-OUT  OF  VARIABLES  ■ 

PRINT (1 -25 X. T. , .F. , 23* . T . >, SUBWGR< . F . > 
PRINTCC1X.TRUE. 

PRINTC  C2)  =  . TRUE . 
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PRINTCC3) = . TRUE .  VAN05050 

PRINTCPPX .TRUE.  VAN05060 

c -  -  -  VAN0507  0 

C -  GROUP  33  MONITOR  PRINT-OUT  i  VAN05080 

C  IXMON<l>,IYMON<l>,IZMON<l>,NPRMON<l>,NPRMNT<l>  VAN05090 

NPRMON=10  VAN05100 

IYM0N=2  VAN05110 

IZM0N=12  VAN05120 

c -  -  -  VAN05130 

C -  GROUP  34  FIELD  PRINT-OUT  CONTROL  :  VAN05140 

C  NPRINT<100>,NTPRIN<100>,NXPRIN<1>,NYPRIN<1>,NZPRIN<1>,  VAN05150 

C  IZPRF<1>,ISTPRF<1>,IZPRL<10000>,ISTPRL<10000>  •  VAN05160 

C  NUMCLS<10>, KOUTPT  VAN05170 

NPRINT=LSWEEP  VAN05180 

c -  -  -  VAN05190 

C -  GROUP  35  TABLE  CONTROL  i  VAN0520D' 

C  TABLES< . F. >, NTABLE, NTABVR, LINTAB, NPRTAB,  NMON,  VAN05210 

C  ITABCl-8),MTABVRCl-8)  VAN05220 

c -  -  -  VAN05230 

C  GROUP  36-38  ARE  NOT  DOCUMENTED  IN  THE  INSTRUCTION  VAN05240 

C  MANUAL  AND  ARE  INTENDED  FOR  MAINTENANCE  PURPOSES  ONLY  VAN05250 

C -  GROUP  36  DEBUG  PRINT-OUT  SLAB  AND  TIME-STEP  i  VAN05260 

C  IZPR1<1>,IZPR2<1>,ISTPR1<1>,ISTPR2<1>  VAN05270 

C -  -  -  VAN05280 

C -  GROUP  37  DEBUG  SWEEP  AND  SUBROUTINES  :  VAN05290 

C  KEMU, KMAIN, KINDEX, KGEOM, KINPUT , KSODAT , KCOMPF, KSORCE,  VAN05300 

C  KSOLV1 , KS0LV2, KS0LV3, KCOMPP,  KADJST , KFLUX, KSHIFT , KDIF,  VAN05310 

C  KCOMPU,  KCOMPV,  KCOMPW, KCOMPR,  KWALL , KDBRH0<-1>, KDBEXP  , KDBMDT  VAN05320 

C  KDBGEN  VAN05330 

c -  - -  -  VAN05340 

C—  GROUP  38  MONITOR,  TEST,  AND  FLAG  >  VAN05350 

C  MONITR<.F.>,FLAG<.F.> , TEST  <.T.>,KFLAG<1>  VAN05360 

C  END  OF  MAINTENANCE-ONLY  SECTION  VAN05370 

c -  -  -  VAN05380 

C -  GROUP  39  ERROR  AND  RESIDUAL  PRINT-OUT  :  VAN05390 

C  IERRP<1000>,RESREFC1, 3-24X25*1.  >,  RESMAP< .  F .  >,  VAN05400 

C  RESI  DC  1-25X2* .  F . ,  23* .  T .  >,  KOUTPT  VAN05410 

RESREFC 1 ) =WIN*RHE  VAN05420 

RESREF(7)=WIN*RESREF(1)  VAN05430 

RESREFC5) =WIN*RESREFC 1 )*0 . 1  VAN05440 

RESREFC HI ) =HTOT*RESREFC 1 )  VAN05450 

RESREFC KE)=RKEIN*RESREFC1)  VAN05460 

RESREFC  EPXEPINXRESREFC1)  VAN05470 

IERRP=LSWEEP/20  VAN05480 

KOUTPT =L SWEEP/ 20  VAN05490 

c -  -  -  VAN05500 

C—  GROUP  40  SPECIAL  DATA  :  LOGICC 1 .  .  10  )  ,  INTGRC 1.  .  10  ) ,  REC  21 .  .30),  VAN05510 

C  NLSP<1>, NISP<1 >, NRSP<1>, SPDATA< .F.>,LSPDAC1),ISPDAC1), RSPDAC 1 )  VAN05520 

C  USE  FIRST  10  ELEMENTS  OF  ARRAYS  LOGIC  8  INTGR  AND  21ST  VAN05530 

C  TO  30TH  OF  ARRAY  RE  FOR  TRANSFERRING  SPECIAL  DATA  FROM  VAN05540 

C  SATELLITE  TO  GROUND,  BUT  IF  REQUIREMENTS  EXCEED  THIS  VAN05550 

C  PROVISION  SET  SPDATA  =  .T.,  AND  DIMENSION  ARRAYS  LSPDA,  VAN05560 

C  ISPDA,  RSPDA  ABOVE  AND  IN  GROUND  AS  NEEDED,  AND  SET  HERE  VAN05570 

c -  VAN05580 

C—  GROUP  42  RESTARTS  AND  DUMPS  i  SAVEM< .  F .  >,  RESTRT< .  F .  >,  KINPUT  VAN05590 

SAVEM= .TRUE.  VAN05600 

BFPLOT a . TRUE.  *  VAN05610 

RESTRT= .TRUE.  VAN05620 

c -  -  -  VAN05630 

c -  -  -  VAN05640 

C -  GROUP  43  GRAFFIC  i  VAN05650 

C  GRAPHS< . F . >, ORTHOG< . T . >, ANTSYM, NPRT<1 >, ITITL<5*4H****>  VAN05660 

C -  FOR  A  GRAFFIC  RUN,  DIMENSION  PHI1  &  PHI2  AS  FOLLOWSi  VAN05670 

C  PHI1CNX*NYXNZ*NM)  VAN05680 

C  PHI2C  C  NX+2) *C NY+2)*C NZ+2)*C  NM+IBLK) )  ,  WHERE  VAN05690 

C  NM=NO.  OF  VARIABLES  STORED  +  DENSITYC-IES)  VAN05700 

C  IBLK=0  IF  BLOCK= . FALSE . , =4  IF  A  3D  RUN,  VAN05710 

C  =3  IF  A  2D.YZ  RUN.  VAN05720 

c -  -  -  VAN057  30’ 

IFCIRUN.EQ.l)  GO  TO  900  VAN05740 

900  CONTINUE  VAN05750 

C -  ALL  RUNS  VAN05760 
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STANDARiTsECTION^STARTSi 
WRITE-GEHERAL~DATA^0N~TOrTHEcGJJ|IE1 ; DTA^ T ISPDA^NISPi RSPDA. HRSP ) 

IFjBS^^>HSftLFC^5:SS^!x^xG^E?Yg^aTz^'IPLAME> 

SND.NX+1,NY+1,NZ+1.NZ,PRTBFC) 

OLD  PRACTICES  RETAINED  FOR  REFERENCE. 

IF(SPDATA)  CALL  SPCDAT(IRUN) 

IF( BLOCK)  CALL  PORDAT(IRUN) 

IF(GRAPHS)  CALL  SORT(IRUN) 

IF(RESTRT)  GO  TO  902 

D°IF(IFIX(FIINIT(INDVAR)+0.1) .NE. 10101)  GO  TO  901 
CALL  FLDDAT ( IRUN) 

GO  TO  902 

901  CONTINUE  ^ 

902  CALL  DATAIO(WRT> 10)  . 

IF(MONITR)  CALL  DATAIO(WRT> -6 ) 

999  CONTINUE 
STOP 
END 

Cxxx  IGEN=1  SO  BFCXYZ  NOT  REQUIRED. 

CXXX  COMMENT  OUT  BOTH  VERSIONS. 


SUBROUTINE  BFCXYZ  ( NXP1 , NYP1 , NZP1 ) 

RETURN 

END 


VAN05770 

VAN05780 

VAN05790 

VAN05800 

VAN05810 

VAN05820 

VAN05830 

VAN05840 

VAN05850 

VAN05860 

VAN0587  0 

VAN05880 

VAN05890 

VAN05900 

VAN059 1 0 

VAN05920 

VAN05930 

VAN05940 

VAN05950 

VAN0596  0 

VAN0597  0 

VAN05980 

VAN05990 

VAN06000 

VAN06010 

VAN06020 

VAN06030 

VAN06040 

VAN06050 
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C$DIRECTIVEx*SATLIT  AMI  LEITNER 

C  LAMINAR  SOLUTION  FOR  NWC5  NY=18  NZ=29  YN=GTH 

C  LECSAT  CONVERTED  TO  DIAMSAT 

C  XFILE  NAME*  MODBFCST . FTN 

C  XABSTRACTi  SATELLITE  MODEL  MAIN  PROGRAM.  THIS  VERSION  IS 
C  FOR  USE  WITH  THE  BODY-FITTED  COORDINATE  SCHEME  (SUMMER  1984 
C  VERSION)  PROVIDED  AS  AN  ATTACHMENT  TO  SPRING  1983  PHOENICS. 

C  ^DOCUMENTATION i  PHOENICS  INSTRUCTION  MANUAL  (SPRING  1983) 

C  WITH  BODY-FITTED  COORDINATES  INSTRUCTION  SUPPLEMENT 

C  (SUMMER  1984). 

C  ^AUXILIARY  SUBROUTINES  (TAPES,  ETC.)  ARE  IN  SATELLITE  LIBRARY 

C  SERVICEU,  WHICH  MUST  BE  INCLUDED  IN  LINK  EDIT  TO  RUN. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  1  STARTS i 


C - 

CHAPTER  1  COMMON  BLOCKS  AND  USER'S  DATA. 

C - - - - - 

INCLUDE  (CMNGUS) 

INCLUDE  (CMNGRF) 

INCLUDE  (GUSSEQ ) 

COMMON/CPI/ IP WRIT , IDUM( 243) 

DIMENSION  GDTAPE(3),DFAULT(4) 

DIMENSION  ARRAY1 (309 ) , ARRAY2( 194) , ARRAY3( 421 ) 

LOGICAL  ARRAY1 , LSPDA, WRT , RD, NAMLST 
INTEGER  ARRAY2,XPLANE, YP LANE, ZP LANE 

INTEGER  P1,PP,U1,U2,V1,V2,W1,W2,R1,R2,RS,EP,H1,H2,H3,C1,C2, 

&real4north,low 

LOGICAL  BFC 

EQUIVALENCE  ( ARRAY1 ( 1 ) , CARTES ) , ( ARRAY2( 1 ) , NX) 

EQUIVALENCE  ( ARRAY3( 1 ) , SPARE1 ( 1 ) ) , ( Ml , R1 ) , (M2, R2) 

EQUIVALENCE  ( LSTRUN, INTGR( 12) ) , (NAMLST , L0GIC(88 ) ) 

EQUIVALENCE  ( LOGIC( 20 ) , BFC) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  1  ENDS. 
C$DIRECTIVE**CMNBF1$$ 

C  THIS  FILE  CONTAINS  SATELLITE  COMMON  BLOCKS  FOR  BFC'S 
C  FI  MUST  BE  DIMENSIONED  TO  GREATER  THAN  OR  EQUAL  TO 

C  (NX+NY+17*NZ+24*NX*NY+6X(NX+1)X(NY+1)+6XND) ’  THE  VALUE 

C  OF  THE  DIMENSION  MUST  BE  SET  AS  NBFC  IN  GROUP  6  OF  SATLIT. 
COMMON/ FOB/ FI ( 5000 ) 

CQMMON/CIB/ND/CIC/KOORD 

COMMON/CID/KDBGG, KDBGMF, KDBGCD, KDBIND, KDBMFX, KDBCDT , KDBPCS, 

8  KDBGUV, KDBGPV 

COMMON/CIE/KDBGS , KDBINS 
COMMON/CI F/ IGEN/CIG/NCART 

C  THE  FOLLOWING  ARRAYS  MUST  BE  EXACTLY  DIMENSIONED  FOR  NXP1, 

C  NYP1  AND  NZP1 ,  BUT  MAY  BE  OVER  DIMENSIONED  FOR  ND. 

C  THE  BFRAC  ARRAYS  MUST  BE  DIMENSIONED  TO  ALLOW  FOR  SETTINGS 
C  IN  SATLIT,  THEY  MAY  BE  OVER  DIMENSIONED. 

COMMON/CRA/XW( 19, 30, 1 )/CRB/XE( 19 , 30,1) 

8  /CRC/YS( 2, 30, 1 )/CRD/YN( 2, 30,1) 

8  /CRE/ZL(2,19,1)/CRF/ZH(2,19,1) 

8  /CRG/RC0N/CRH/DARCY/CRI/BXFRAC(99)/CRJ/BYFRAC(99) 

8  /CRK/BZFRAC(99) 

COMMON/CL A/ST0RSA(6 ) , ST0RWD(6 ) , STORP, STORPE, STORPN, 

8  STORPH, ST0R1 , ST0R2, ST0R3 , STOUNV, PRTBFC, STOCRN 

COMMON/CLC/BFPLOT 

LOGICAL  STORP, STORPE, STORPN, STORPH, ST0R1 , ST0R2, ST0R3, 

8  STORSA,STORWD, STOUNV, PRTBFC, BFPLOT, STOCRN 

C  END 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  USER  SECTION  1  STARTS: 

C  GRAFFIC  ARRAYS  DIMENSIONED  AS  NEEDED... 

C0MM0N/GRAF1/PHIK1)  /GRAF2/PHI2(  1 ) 

C  POROSITY  8  SPECIAL  DATA  ARRAYS  DIMENSIONED  AS  NEEDED... 
DIMENSION  PE(1,1,1),PN(1,1,1),PH(1,1,1),PC( 1,1,1) 

DIMENSION  LSPDA(1),ISPDA(1),RSPDA(1) 

C  USER  PLACES  HIS  VARIABLES,  ARRAYS,  EQUIVALENCES  ETC.  HERE. 
EQUIVALENCE( RAIR, RE( 21 ) ) , ( GAMA, RE( 22) ) , ( GSWP , RE( 23 ) ) 

1 , (GPR, RE( 24 ) ) , (TW, RE( 25 ) ) , (GEMU1, RE(26 ) ) , ( JEMU1 , INTGRC 1 ) ) 

C  USER  PLACES  HIS  DATA  STATEMENTS  HERE. 

DATA  NLSP,NISP, NRSP/1, 1,1/ 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  USER  SECTION  1  ENDS. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  2  STARTS: 


VAN00010 
VAN00020 
VAN00030 
VAN00040 
VAN00050 
VAN0006  0 
VAN0007  0 
VAN00080 
VAN00090 
VAN00100 
VAN001 10 
VAN00120 
VAN00130 
VAN00140 
VAN00150 
VAN0016  0 
VAN00170 
VAN00180 
VAN00190 
VAN00200 
VAN0021 0 
VAN00220 
VAN00230 
VAN00240 
VAN00250 
VAN0026  0 
VAN0027  0 
VAN00280 
VAN00290 
VAN00300 
VAN00310 
VAN00320 
VAN00330 
VAN00340 
VAN00350 
VAN0036  0 
VAN00370 
VAN00380 
VAN00390 
VAN00400 
VAN0041 0 
VAN00420 
VAN00430 
VAN00440 
VAN00450 
VAN0046  0 
VAN0047  0 
VAN00480 
VAN00490 
VAN00500 
VAN00510 
VAN00520 
VAN00530 
VAN00540 
VAN00550 
VAN0056  0 
VAN0057  0 
VAN00580 
VAN00590 
VAN00600 
VAN006 1 0 
VAN00620 
VAN00630 
VAN00640 
VAN00650 
VAN0066  0 
VAN0067  0 
VAN00680 
VAN00690 
VAN007  00 
VAN007 10 
VAN00720 
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CHAPTER  2  SET  CONSTANTS,  AND  ARRANGE  FILE  MANIPULATIONS. 


C- 

C 

C 


CD 

CD 


PLEASE  DO  NOT  ALTER,  OR  RE-SET,  ANY  OF  THE  REMAINING 
STATEMENTS  OF  THIS  CHAPTER, 

DATA  CELL, EAST, WEST, NORTH, SOUTH, HIGH, LOW, VOLUME/ 

&  0.,1.,2.,3.,4.,5.,6.,7 .  / 

DATA  P1,PP,U1,U2,V1,V2,W1,W2,R1,R2,RS,KE,EP,H1,H2,H3,C1,C2, 
8C3,C4/1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20/ 

DATA  FIXFLU,FIXVAL,ONLYMS,WALL/1.E-10,1 .E10, 0.0, -10.0/ 

DATA  IPLANE,XPLANE,YPLANE,ZPLANE/0, 1,2,3/ 

DATA  WRT, RD,DFAULT/. TRUE. , .FALSE. , 4HDEFA, 4HULT . /4HDTA/ , 1HG/ 
DATA  GDTAPE/4HGUSI ,4HE1.D,2HTA/ 

DATA  NL DATA, NI DATA, NRDATA/3 09, 194, 421/ 

DATA  NLCREG, NTCVRG/6  0,350/ 

DATA  TITPP,TITC1,TITC2,TITC3/3HRH0, 4HMACH , 4HTEMP , 4HCFST/ 

CALL  TAPES(10,GDTAPE,3,1,4*NRDATA) 

- READ  DEFAULT  FILE  IF  BLOCKDATA  ABSENT 

IF(INTGR1(29) .NE. 10)  GO  TO  2 

CALL  WRIT40C  40HDATA  ESTABLISHED  IN  BLOCK  DATA.  ) 

GO  TO  3 
2  CALL  DEFLT 

2  CALL  TAPES(1,DFAULT,4,2,4*NRDATA) 

CALL  DATAIOCRD, 1 ) 

CALL  WRIT40C  40HDATA  TAKEN  FROM  DEFAULT. DTA  ON  GROUP  A/C) 

3  CALL  WRIT40C40HFILE  MODSTL.FTN  IS  THE  SATLIT  USED.  ) 

LOGICC89 ) = . TRUE . 


C - 

CHAPTER  3  DEFINE  DATA  FOR  NRUN  RUNS. 

C - 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  2  ENDS. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  USER  SECTION  2  STARTS* 

C—  GROUP  41MULTI-RUNS  *  RUN(  1-30  )<  .T  . ,  29* .  F .  > 

C 

RUN(1)=. FALSE. 

C  NOTE:  ALL  RUNS  ARE  DEACTIVATED  AT  THIS  POINT  -  USER  SHOULD 
C  ====  SWITCH  ON  ONE  ONLY  OF  RUNS  1-4  IN  NEXT  STATEMENT. 

RUN( 1 ) = . TRUE . 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  USER  SECTION  2  ENDS. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  3  STARTS* 
DO  10  IRUN=1,30 

I F C . NOT . RUN( I  RUN  ) )  GO  TO  10 
NRUN=NRUN+1 
LSTRUN=IRUN 
10  CONTINUE 

DO  999  IRUN” 1 , L STRUN 

IF( .NOT  RUN(IRUN))  GO  TO  999 
INTGR(ll)  =  IRUN 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  3  ENDS. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  USER  SECTION  3  STARTS* 

C—  ALL  INTEGER  VARIABLES  ARE  DEFAULTED  TO  0,  AND  REAL  VARIABLES 
C  TO  0.0,  UNLESS  OTHERWISE  INDICATED. 

C  E.G.  BY  VARIABL E<1 0>,  OR  <10. 0>  AS  APPROPRIATE. 

C  THE  DEFAULT  SETTINGS  OF  ALL  LOGICAL  VARIABLES  ARE  ALWAYS 
C  INDICATED,  E.G.  VARIABLE< . T . >,  OR  VARIABL E< . F . > . 

C 

C -  RUN1 

C . . —  -  - 

C -  GROUP  1.  FLOW  TYPE  * 

C  PARAB< . F . >, CARTES< . T . >, ONEPHS< . T . > 

C -  -  - 

C—  GROUP  2.  TRANSIENCE  * 

C  .  STEADY< . T . >, ATIME, LSTEP<1>, FSTEP<1> 

C  TLAST<1  .El  0>,TFRAC(  1-30X30X1  .> 

C  SERVICE  SUBROUTINE  FOR  f NT f  POWER-LAW  TIME  STEPS* 

C  CALL  GRDPWRCO, NT, TLAST, POWER) 

C -  -  - 

C -  GROUP  3.  X-DIRECTION  * 

C  NX<1>,XULAST<1 .0>, XFRACC 1-30 ) 

C  SERVICE  SUBROUTINE  FOR  POWER-LAW  GRID* 

C  CALL  GRDPWRC1, NX, XULAST, POWER) 

C -  -  - 


VAN00730 
VAN007  40 
VAN007  50 
VAN00760 
VAN0077  0 
VAN00780 
VAN007  90 
VAN008  00 
VAN008 1 0 
VAN00820 
VAN00830 
VAN00840 
VAN00850 
VAN0086  0 
VAN0087  0 
VAN00880 
VAN00890 
VANO  09  00 
VAN009 1 0 
VAN00920 
VAN00930 
VANO  0940 
VAN00950 
VAN0096  0 
VANO  097  0 
VAN0098  0 
VANO  099  0 
VANO 1 000 
VAN01010 
VANO 1020 
VANO 1 030 
VANO 1 040 
VANO 1050 
VAN01060 
VANO 1 07  0 
VAN01080 
VAN01090 
VANO 1 100 
VAN01110 
VANO 1120 
VANO 1130 
VANO 1 140 
VANO 1 1 50 
VANO 1 16  0 
VAN01170 
VANO 1 18  0 
VAN01190 
VANO 1 200 
VANO 1 21 0 
VAN 01220 
VAN01230 
VAN 01240 
VANO 1250 
VAN0126  0 
VAN0127  0 
VAN01280 
VAN 01290 
VANO 1 30  0 
VAN01310 
VANO 1 320 
VAN01330 
VANO 1 340 
VAN01350 
VAN0136  0 
VANO 1 37  0 
VANO 1380 
VAN01390 
VAN01400 
VAN01410 
VANO 1 420 
VANO 1 43  0 
VANO 1 440 
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—  GROUP  4.  Y-DIRECTION  i 

NY<1>,YVLAST<1.0>,YFRAC(1-30),RINNER,SNALFA 
SERVICE  SUBROUTINE  FOR  POWER-LAW  GRID« 

CALL  GRDPWR(2, NY, YVLAST, POWER) 

NY=18 


—  GROUP  5.  Z-DIRECTION  i 

NZ<l>,ZWLAST<1.0>,ZFRAC(l-30) 

SERVICE  SUBROUTINE  FOR  POWER-LAW  GRID: 
CALL  GRDPWR(3,NZ,ZWLAST, POWER) 

NZ=29 


—  GROUP  6.  MOVING  GRID  OR  DISTORTED  (BODY-FITTED)  GRID  : 

—  MOVING  GRID  : 

MGRID, IZW1 , IZW2, AZW2, BZW2,CZW2, PINT,ZW2M1T 

-  BODY-FITTED  GRID  — 

BFC< .  T .  >, IGEN<1>,ND<1>,NBFC<5000>,KQQRD, RCON 
BXFRACC  l-NXXl .  0  ,  NXM1X0  .  0> 

BYFRACC  l-NYXl .  0 ,  NYMlxO  .  0> 

BZFRACC 1-NZ) <1 . 0/ NZM1X0 . 0> 

SERVICE  SUBROUTINE  FOR  SUB-DOMAIN  SPECIFICATION  (FOR  IGEN=1 
ONLY) * 

CALL  DOMAIN(ID,IXF, IXL , IYF, IYL , IZF, IZL) 
XE(l-NYPl,l-NZPl,l-NDX(NYPlxNZPlxND)xl .0>, 

XW( 1-NYP1 , 1-NZP1 , 1 -ND) , 

YN(1-NXP1,1-NZP1,1-ND)<(NXP1XNZP1*ND)X1.0>, 

YS( 1-NXP1 , 1-NZP1 , 1-ND) , 

ZH(l-NXPl,l-NYPl,l-NDX(NXPlxNYPlxND)xl .0>, 

ZL(1-NXP1, 1-NYP1,1-ND), STORSA( 1-6X6 X . F . > , STORWDC 1-6) <6x.F.>, 
STORP<.F.>,STORPE<.F.>,STORPN<.F.>,STQRPH<.F.>,STOUNV<.F.>, 
PRTBFC< . F . >,  DARCY , BFPLOT< . F . > 

CYCLIC  BOUNDARY  CONDITIONS  ARE  DEFAULTED  INACTIVE  ; 

TO  ACTIVATE  THEM  AT  SELECTED  IZ  SLABS  USE  SERVICE  SUBROUTINE* 
CALL  XCYIZCIZ, .TRUE. ) 

SERVICE  SUBROUTINE  TO  DEACTIVATE  CURVATURE  TERMS  IN  U,  V 
AND  W  EQUATIONS  ASSOCIATED  WITH  CURVATURE  OF  IX,  IY,  IZ 
GRID  LINES  RESPECTIVELY* 

CALL  UCURVECIZ, .FALSE. ) 

CALL  VCURVE(IZ, .FALSE.) 

CALL  WCURVE(IZ, .FALSE. ) 

NCART<1> 

XWARNINGSI I | | | | 

A)  WHEN  USING  BFCS  ST0VARCH3) ,  STOVARC C4) ,  ST0VARC21 )  ARE 
AVAILABLE  ONLY  FOR  STORING  NON-ORTHOGONAL  VELOCITY 
COMPONENTS. 

B)  MULTI-RUNS  ARE  NOT  ALLOWED  WITH  BFC  OPTION. 

C)  MOVING  GRID, TWO-PHASE  AND  PARABOLIC  OPTIONS  ARE  NOT 
AVAILABLE  WITH  BFC  OPTION. 

D)  KE-EP  TURBULENCE  MODEL  SHOULD  BE  USED  WITH  BFC'S  ONLY 
WHEN  THE  MAIN  FLOW  IS  IN  THE  IZ  DIRECTION. 

E)  BUILT-IN  GRAVITY  TERMS  DO  NOT  TAKE  ACCOUNT  OF  BFC’S. 
XNOTES 


A)  THE  STANDARD  VELOCITY-FIELD  PRINTOUT  FOR  THE 
VELOCITY  RESOLUTES  IS  ACTIVATED  JN  THE  USUAL 
WAY.  AN  ADDITIONAL  OPTION  EXISTS  FOR  PRINTING  THE 
CARTESIAN  VELOCITY-COMPONENTS  WHICH  MAY  BE 
ACTIVATED  BY  SETTING  THE  FOLLOWING  LOGICALS: 

STOVARC U2)  =  . T .  FOR  U-CQMPQNENT  (CARTESIAN) 
ST0VAR(V2X.T.  FOR  V-COMPONENT  (CARTESIAN) 
ST0VAR(W2)=.T.  FOR  W-COMPQNENT  (CARTESIAN) 

SIMILARLY  PRINTOUT  OF  NON-ORTHOGONAL  VELOCITY 
COMPONENTS  MAY  BE  ACTIVATED  AS  FOLLOWS* 

STQVAR(C4 ) = . T .  FOR  U-COMPQNENT  (NQN-QRTHOG) 
ST0VAR(H3)=.T.  FOR  V-COMPONENT  (NON-ORTHOG) 

STQVAR( 21 ) = . T .  FOR  W-COMPQNENT  (NON-ORTHOG) 

B)  BFC  (TO  ACTIVATE  THE  BFC  OPTION),  IGEN  (THE  CODE  FOR  METHOD 

OF  GRID  SPECIFICATION),  ND  (NUMBER  OF  SUB-DOMAINS)  AND 
NBFC  (THE  FI  ARRAY  DIMENSION),  MUST  BE  SET  BEFORE 
"STANDARD  BFC  SECTION  2".  ====== 


VAN01450 

VAN0146  0 

VAN01470 

VANO 1480 

VAN01490 

VAN01500 

VAN01510 

VAN01520 

VAN01530 

VAN01540 

VAN01550 

VAN01560 

VAN0157  0 

VAN01580 

VAN01590 

VAN016DO 

VAN016 1 0 

VANO 1620 

VAN01630 

VAN01640 

VAN01650 

VAN0166  0 

VAN01670 

VAN01680 

VAN01690 

VAN017  00 

VAN017 10 

VAN01720 

VAN017  30 

VAN01740 

VAN01750 

VANO 176  0 

VAN0177  0 

VAN01780 

VAN017  90 

VAN01800 

VAN01810 

VAN01820 

VAN01830 

VAN01840 

VAN01850 

VANO 186  0 

VAN01870 

VANO 1880 

VAN01890 

VAN01900 

VAN0191 0 

VAN01920 

VAN01930 

VAN01940 

VAN01950 

VAN01960 

VAN01970 

VAN01980 

VAN01990 

VAN02000 

VAN0201 0 

VAN02020 

VAN020  30 

VAN02040 

VAN02050 

VAN0206  0 

VAN0207  0 

VAN02080 

VAN02090 

VAN02100 

VAN021 1 0 

VAN02120 

VAN021 30 

VAN02140 

VAN021 50 

VAN0216  0 


50 


FILE*  VAN4SAT  FORTRAN  A1 


C  ALL  OTHER  BFC  DATA  MUST  BE  SET  AFTER  "STANDARD  BFC 

C  SECTION  2.  ===== 

r  C)  NXP1,  NYP1,  NZP1  STORE  NX+1,  NY+1,  NZ+1;  THESE  ARE 

C  AVAILABLE  TO  USER  AFTER  STANDARD  BFC  SECTION  2. 

C  D)  FOR  IGEN=1  USE  BXFRAC, BYFRAC  &  BZFRAC  IN  PLACE  OF 

C  XFRAC, YFRAC  8  ZFRAC . 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  BFC  SECTION  1  STARTS: 
C  DEFAULT  SETTINGS* 

NCART=10 
BFC= . TRUE. 

IGEN=1 
ND  =  1 

NBFC=5000 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  BFC  SECTION  1  ENDS.: 

C  XUSER  SETS  BFC#  IGEN#  ND  AND  NBFC  HERE* 

C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  BFC  SECTION  2  STARTS* 
CALL  SB4I CNXP1, NX+1, NYP1, NY+1, NZP1, NZ+1, I#0) 

IF(BFC)  CALL  BFCDFTC NBFC, XE, XW, YN# YS,ZH, ZL , ND, NXP1 , NYP1 , 

&  NZP1 , NZ) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  BFC  SECTION  2  ENDS. 

C  XUSER  SETS  ALL  OTHER  BFC  VARIABLES  HERE* 

C  XUSING  NONIFORM  GRID  1-8 
GTH=65 . E-3 
GTL  =  150 . E-3 
GBETA=4 . 

GBETA=GBETAX3. 1415927/180 
GTAB=TAN(GBETA) 

DELMAX=2.E-3 
GNBL  =  5 . 

GPWR=4 . 

DO  64  IY=1 , 5 

64  BYFRAC(IY)=(FLOAT(IY)/GNBL)XXGPWRXDELMAX/GTH 
BYFRAC (6 )= BYFRAC C 5 )+3 . E-3/GTH 

DEL=( 1 . - BYFRAC ( 6 ) )/ ( FLOAT ( NY) -GNBL -1 ) 

DO  65  IY=7 , NY 

65  BYFRAC ( IY)=BYFRAC( IY-1 )+DEL 

BZFRAC ( 1 ) =1 0 . E-3 
DO  66  IZ=2, 5 

66  BZFRACC IZ) =10 . E-3+BZFRACC IZ-1) 

BZFRACC  6 ) =BZFRACC  5)  +  5 . E-3 

DO  67  IZ=7 , 9 

67  BZFRACC IZ)=BZFRACC IZ-1 )+2. E-3 
DO  68  IZ=10,11 

68  BZFRACC IZ) =BZP RAC C IZ-1 ) + . 5E-3 
BZFRACC 12) =BZFRACC 11 )+l .E-3 
DO  77  IZ=13, 14 

77  BZFRACC IZ)=BZFRACC IZ-1 )+.5E-3 
DO  78  IZ=15, 15 

78  BZFRACC IZ)=BZFRACC IZ-1 )+l. E-3 
BZFRACC 16 )=BZFRACC1 5 )+l .E-3 
BZFRACC 17 ) =BZFRACC 16 )+2 . E-3 
BZFRACC 18 )=BZFRACC 17 )+7. E-3 
DO  69  I Z= 19 , 22 

69  BZFRACC IZ)=BZFRACC IZ-1 )+10. E-3 
BZFRACC 23)=BZFRACC 22 )+3. E-3 
BZFRACC  24 )  =  BZFRACC  23 ) +  2 . E-3 
BZFRACC 25) =BZFRACC 24 )+2. E-3 
BZFRACC  26 )  =  BZFRACC  25 )  +  3 . E-3 
BZFRACC  27 ) =BZFRACC  26 )  +  5 . E-3 

DO  71  IZ=28 , NZ 

71  BZFRACC I Z)=BZFRACC IZ-1 )+10. E-3 
DO  72  IZ=1 , NZ 

72  BZFRACC IZ)=BZFRACCIZ)/GTL 
CALL  DOMAIN C 1,1,NX,1,NY,1#NZ) 

DO  61  IX=1 , NXP1 

DO  62  I Y=1 , NYP1 
ZL  C IX, IY, 1 )  =  0 . 0 
62  ZHC IX, IY, 1 ) =GTL 
DO  63  IZ=1 , NZP1 
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YNCIX,IZ,1)=GTH 
63  YSCIX,IZ,1)=0.0 

C  YSCIX,13,1)  SHOULD  COME  AFTER 
DO  662  IZ=5, 25 
CBL  DO  662  IZ=16,25 

662  YS(IX,IZ,l)=(BZFRAC(IZ-l)-BZFRAC(3))xGTAB*GTL 
CBL  DO  663  IZ=13,15 

CBL  GZ12  =  C  BZFRACC IZ-1 ) -BZFRACC 1 1 ) )XGTL- . 5E-3 

CBL663YS(IX,IZ,1)=SQRT(YSCIX,16,1)*GZ12*2.-GZ12*X2) 

DO  664  IZ=26,NZ 
664  YS(IX,IZ,1)=YSCIX,25,1) 

61  CONTINUE 

STORSACIFIXC LOW) )=. TRUE. 

ST ORSA(IFIX( HIGH))3. TRUE. 

STORSACIFIXC SOUTH))3. TRUE.  . 

STORWDC I FIXC SOUTH) ) 3 . TRUE.  ' 

STORP3 . TRUE . 

PRTBFC=.TRUE. 

CDAR  DARCY=1 . E10 

C -  -  - 

C—  GROUP  7.  BLOCKAGE)  BLOCK< . F .>, IPLANE, IPWRIT 
C  XSET  CONSTANT  POROSITIES  OVER  SUB-DOMAINS  USING) 

C  CALL  CONPORCIR, TYPE, VALUE, IXF, IXL , IYF, IYL , IZF, IZL ) ,  WHERE) 

C  IR=RUN  SECTION  NUMBER,  E.G.  1  FOR  RUN1  SECTION;  'TYPE'3  EAST, 

C  WEST,  NORTH,  SOUTH,  HIGH,  LOW  &  CELL.  'VALUE' =WANTED  POROSITY 

C  OVER  REGION  IXF, . . . IZL . 

C  xDIMENSION  ARRAYS  PECNX, NY, NZ) ,  PNCNX, NY, NZ) ,  PHCNX, NY, NZ) ,  & 

C  PCCNX,NY,NZ)  ABOVE. 

C  XFOR  FULLY-BLOCKED  CELLS  CIE.  'VALUE'3  0.0)  USER  NEED  SET  ONLY 

C  THE  'CELL'  POROSITY  CTO  ZERO),  AS  CELL-FACE  AREAS  ARE  THEN 

C  AUTOMATICALLY  ZEROED. 

C  XFOR  SATELLITE  PRINTOUT  OF  ALL  POROSITIES  IN  DOMAIN,  'IPLANE'3 
C  XPLANE  YPLANE  OR  ZPLANE,  FOR  DESIRED  CROSS-SECTION  DIRECTION. 
C  XFOR  EACH  'TYPE'  A  MAXIMUM  OF  10  CALLS  TO  CONPOR  IS  ALLOWED, 

C  BUT  IF  REQUIREMENTS  EXCEED  THIS  PROVISION  SET  BLOCK3. T.  & 

C  IPWRIT  =  -l ,  AND  SET  POROSITY  ARRAYS  EXPLICITLY  HERE  AS  WANTED. 

C  IN  THIS  CASE,  THE  USER  MUST  SET  A  L  L  ELEMENTS  OF 

C  ARRAYS  PE,  PN,  PH,  PC  CMANY  MAY  BE  0.0  OR  1.0).  HE  MAY  USE) 

C  CALL  CRCPARRAY, VALUE, IXF, IXL, IYF, IYL, IZF, IZL, NX, NY, NZ) 

C  ANY  NUMBER  OF  TIMES,  TO  SET  'PARRAY'  C=  PE,  ETC.)  TO 
C  'VALUE'  OVER  RANGE  IXF  TO  IXL,  IYF  TO  IYL,  IZF  TO  IZL. 

C  XCONPOR  MUST  NOT  BE  USED  IN  CONJUNCTION  WITH  EXPLICIT 
C  SETTINGS  OF  THE  ARRAYS  CINCLUDING  SETTINGS  VIA  CR) . 

C -  -  - 

C -  GROUP  8. DEPENDENT  VARIABLES  TO  BE  SOLVED  FOR  OR  STORED  > 

C  SOL  VARC 1-25  )  <25X.  F.  >,  STOVARC 1-25X25X .  F .  >,  CONC1 C  1-4  )<4X.T.> 

C  USE  FOLLOWING  NAMED  INTEGERS  FOR  ARRAY  ELEMENTS  1-20) 

C  P1,PP,U1,U2,V1,V2,W1,W2,M1,M2,RS,KE,EP,H1,H2,H3,C1,C2,C3,C4. 

S0LVARCP1)=.TRUE. 

SOLVARCPP)=.TRUE. 

S0LVARCV1)=.TRUE. 

S0LVARCW1)=.TRUE. 

S0LVARCH1)=.TRUE. 

CT  SOLVARCKE)=.TRUE. 

CT  SOLVARCEP)=.TRUE. 

ST0VARCV2)=.TRUE. 

ST0VARCW2)=.TRUE. 

ST0VARCC1 )= .TRUE . 

ST0VARCC2) 3 . TRUE . 

ST0VARCC3)3 .TRUE. 

C -  -  - 

C -  GROUP  9.  VARIABLE  LABELS  > 

C  TITLEC1-25X2HP1,2HPP,2HU1,2HU2,  2HV1 ,2HV2,2HW1,2HW2,2HR1, 

C  2HR2, 2HRS, 2HKE, 2HEP, 2HH1 , 2HH2, 2HH3, 2HC1 , 2HC2, 

C  2HC3 , 2HC4 , 2HRX, 2HRY , 2HRZ,  2x4Hxxxx> 

TITLECC1)=TITC1 

TITLECC2)=TITC2 

TITLECC3)=TITC3 

TITLECPP)=TITPP 

C -  -  - 

C—  GROUP  10  PROPERTIES) 

C  IRH01<1>,IRH02<1>,RH01<1 .0>,RH02<1 .0>, 
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ARH01<1 .0>,BRH01<1.Q>,CRH01<1 .0> 
IEMU1<1>,EMU1<1.0>,EMULAM<1.E-10> 
IHSAT,H1SAT,H2SAT,PSATEX<1 .0> 

SIGMA<1-25X1.0,2.0,1.,1.E10,1.#1.E10,1.,1.E10, 
4X1. 0,1. 314, 1.0,1. El  0, 10X1. 0> 

IRH01=-1 
PT0T=55.E5 
T0T=555 . 55 
RAIR=287. 

GAMA=1 .35 

CP=RAIR/ ( 1-1/GAMA ) 

TW=323 . 

HWALL=TWXCP 

htot=cpxtot 

RHT0T=PT0T/T0T/RAIR 
L0GICC87 ) = . TRUE . 

ARH01=RHT0T/PT0TXX( 1/GAMA) 

BRH01=1 ./GAMA 
TURBULENT  OR  LAMINAR 
I EMU1 =  -l 
I EMU1  = 1 
JEMU1 =1 EMU1 
EMU1 =1 . E-5 
EMULAM=EMU1 
GEMU1 =EMU1 
GPR  = . 7 

SIGMAC  24 ) =GPR 
SIGMAC 14 ) =GPR 


-  GROUP  11  INTER-PHASE  TRANSFER  PROCESSES  * 

ICFIP,CFIPS,IMD0T,CMD0T,CA1I<1 .E6>,CA2I<1 .E6> 


—  GROUP  12  SPECIAL  SOURCES  * 

ISPCSOC 1 -25 ) , AGRAVX, AGRAVY, AGRAVZ, ABUOY , HREF 


—  GROUP  13  INITIAL  FIELDS  * 

FI  I N IT ( 1-25X25X1 .  E-10> 

MACH  NO.  OF  FREE  STREAM 
GMACH=3 . 2 

A=1+(GAMA-1)/2XGMACHXX2 

TE=TOT/A 

RHE=RHTOT/AXX(l/( GAMA-1)) 

PSTAT=PTOT/AXX( GAMA/ ( GAMA -1 ) ) 

RH01=ARH01XP$TATXXBRH01 

SONIC=SQRT(GAMAXRAIRXTE) 

WIN=SONICXGMACH 

RKEIN=0.01XWINXX2 

EPIN=0 . 16XRKEINXX1 . 5/GTH/2 . 

FIINITCW1XWIN 

FIINIT(P1)=PSTAT 

FI INIT ( HI ) =HTOT 

FIINIT(KE)=RKEIN 

FI INIT ( EP ) =EPIN 


-  GROUP  14  BOUNDARY/INTERNAL  CONDITIONS  * 

ILOOP1, ILOOPN,XCYCLE< . F .>, PBAR, REGION! 1-10 XI  OX . T . > 

XN.B.  ALL  10  REGIONS  ARE  DEFAULTED  .TRUE..  THE  USER  SHOULD 
SET  REGION(I)=. FALSE.  FOR  UNUSED  REGIONS  fIf. 

DO  14  1=1,10 
14  REGION(I)=. FALSE. 


-  GROUP  15  TO  24;  REGIONS  1  TO  10 

ONLY  THOSE  REGIONS  ARE  ACTIVE  WHICH  ARE  SPECIFIED  BY  THE 
USER,  PREFERABLY  BY  WAY  0F«- 

CALL  PLACE(IREGN,TYPE,IXF,IXL,IYF,IYL,IZF,IZL)  8 
CALL  COVALC IREGN, VARBLE,COEFF, VALUE) 

CALL  PLACE(l,LOW,l,NX,l,NY,l,l) 

C  CALL  C0VAL(1,M1, FIXFLU, WINXRHE) 

CDAR  CALL  COVAL ( 1 , Ml , 1 . E-20 , 1 . E+20XWI NXRHE) 
GCM=2XGAMA/WIN/(GAMA-1 ) 

GVM=PTOTXRHE/RHTOT 
CALL  C0VAL(1,M1,GCM,GVM) 
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C 

C 

C 


C 


C 

C 


CT 

CT 

C-- 

C— 

C 

C 

C— 

C— 

C 


CALL  COVAL(l,Wl,ONLYMS,WIN) 

CALL  COVAL ( 1  >  HI  > ONLYMS, HTOT) 

CALL  COVAL ( 1 »KE»  ONLYMS, RKEIN) 

CALL  COVAL ( 1 > EP>  ONLYMS, EPIN) 

CALL  PLACE(2,HIGH,1,NX,1,NY,NZ,NZ) 

CALL  COVAL(2,M1,FIXVAL,PSTATXO. ) 

CALL  COVALC2,M1,1000XWINXRHE/PSTAT,PSTAT) 
CALL  C0VAL(2,H1,0NLYMS,HT0T) 

WALL  ALONG  THE  VANE  IZ(11,NZ) 

GCM=EMUl/( ,5XBYFRAC(1)XGTH) 

DY1=BYFRAC( 1 )xGTH 

GOEFF=EMU1/(0.5XDY1) 

GOEFH=EMUl/(0.5xDYlxSIGMA(24)) 

CALL  PLACE(3,S0UTH,1,NX,1,1,4,NZ) 

CALL  COVAL(3,W1,GOEFF,0. ) 

CALL  C0VALC3,H1,G0EFH,HWALL)- 
CALL  COVAL(3,W1,WALL,0. ) 

CALL  COVAL (3, HI , WALL ,  HWALL ) 

CALL  COVALC3,KE,WALL,0. ) 

CALL  COVAL(3,EP,WALL,0. ) 


GROUP  25  GROUND  STATION  . 

GROSTA< . F . >, NAMLST < . F .  > 

XNAMLST  ACTIVATES  NAMELIST  IN  GROUND. 
GROSTA= . TRUE . 


GROUP  26  SOLUTION  TYPE  AND  RELATED  PARAMETERS 
WHOL EP< . F . > , SUBPST < . F . >, DONACC<  .  F.  > 
WHOLEP=.TRUE. 


t 


c -  -  - 

c—  GROUP  27  SWEEP  AND  ITERATION  NUMBERS  » 

C  FSWEEP<1>,LSWEEP<1>,LITHYD<1>,LITC<1>,LITKE<1>,LITH<1>, 
C  L ITERC 1-25 )<9*1, -1,1 5*1> 

C  IVELF<1>,NVEL<1>,IVELL<10000>, 

C  IKEF<1>,NKE<1>,IKEL<10000>, 

C  IENTF<1>,  NENT<1>, IENTL<10000># 

C  ICNCF<1>,NCNC<1>,ICNCL<10000>, 

C  IRH01F<1>,NRH01<1>,IRH01L<10000>, 

C  IRH02F<1>,NRH02<1>,IRH02L<10000> 

LSWEEP=400 
GSWP=LSWEEP 
CR  F$WEEP=200 

LITER(PP) =20 
LITERC VI ) =5 
L ITERC W1 ) =5 
C  LITHYD=2 

C -  -  - 

C—  GROUP  28  TERMINATION  CRITERIA  i 

C  ENDIT  (1-25X9^1  .  E-10 , 0 . 5 , 15*1  .  E-10> 

ENDIT ( 1 ) =1 . E-5 


C -  -  - 

C—  GROUP  29  RELAXATION  : 

C  RLXP<1 . >, RLXPXY<1 . > , RLXPZ<1 . > , RLXRH0<1 . > , RLXMDT<1 . >  > 

C  DTFAL S( 3-25) <23X1. El 0> 

DTFALS(Wl)=l.E-5 
DTFALS(V1)=1. E-5 
RLXP= . 2 

C -  -  - 

C—  GROUP  30  LIMITS  » 

C  VELMAX<1 .E10>,VELMIN<-1 . E10>, RH0MAX<1 . E10> , RH0MIN<1 .E-10>, 
C  TKEMAX<1 .E10>,TKEMIN<1.E-10>, EMUMAX<1 . E10>, EMUMIN<1 . E-10>, 
C  EPSMAX<1 . E10>, EPSMIN<1 . E-l 0>, AMDTMX<1 . El 0> , AMDTMN<-1 . E10> 

C - T -  -  - 

C -  GROUP  31  SLOWING  DEVICES  *  SL0RH0<1 . >, SL0EMU<1 . > 

SLORHO= . 2 

C - -  -  - 

C — -  GROUP  32  PRINT-OUT  OF  VARIABLES  i 
C  PRINT ( 1-25) < . T , , . F. , 23x .T . >, SUBWGR< . F . > 

PRINTCC1 ) = . TRUE . 

PRINT(C2)=.TRUE. 

PRINT(C3)=.TRUE. 

PRINT(PP)=.TRUE. 
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—  GROUP  33  MONITOR  PRINT-OUT  * 

IXM0N<1>,IYM0N<1>,IZM0N<1>,NPRM0N<1>,NPRMNT<1> 

NPRM0N=5 

IYM0N=2 

IZM0N=12 


-  GROUP  34  FIELD  PRINT-OUT  CONTROL  « 

NPRINT<100>,NTPRIN<100>,NXPRIN<1>,NYPRIN<1>,NZPRIN<1>, 
IZPRF<1>,ISTPRF<1>, IZPRL<10000>,ISTPRL<10000> 
NUMCLS<10>, KOUTPT 
NPRINT=LSWEEP 


-  GROUP  35  TABLE  CONTROL  1 

TABLES<.F.>,NTABLE,NTABVR,L INTAB, NPRTAB,NMON, 
ITABC1-8) , MTABVRC 1-8 ) 


GROUP  36-38  ARE  NOT  DOCUMENTED  IN  THE  INSTRUCTION 
MANUAL  AND  ARE  INTENDED  FOR  MAINTENANCE  PURPOSES  ONLY 
—  GROUP  36  DEBUG  PRINT-OUT  SLAB  AND  TIME-STEP  1 
IZPR1<1>, IZPR2<1>, ISTPR1<1>,ISTPR2<1> 


-  GROUP  37  DEBUG  SWEEP  AND  SUBROUTINES  1 

KEMU, KMAIN, KINDEX, KGEOM, KINPUT, KSODAT , KCQMPF, KSQRCE, 
KS0LV1,KS0LV2,KS0LV3, KCOMPP , KADJST,KFLUX,KSHIFT,KDIF, 
KCOMPU, KCOMPV, KCOMPW, KCOMPR, KWALL , KDBRH0<-1>, KDBEXP, KDBMDT 
KDBGEN 


-  GROUP  38  MONITOR, TEST, AND  FLAG  * 

MONITR<.F.>, FLAG<.F.>,TEST<.T. >,KFLAG<1> 
END  OF  MAINTENANCE-ONLY  SECTION 


—  GROUP  39  ERROR  AND  RESIDUAL  PRINT-OUT  t 

I ERRP<1000>,RESREF( 1,3-24) <25*1. >,RESMAP<.F.>, 
RESI DC  1-25 ) <2*. F. ,23X.T.>, KOUTPT 
RESREFC 1 ) =WINXRHE 
RESREF(7)=WINXRESREF(1) 

RESREFC 5 )=WINXRESREFC1)X0.1 
RESREFC HI )=HT0TXRESREF(1) 

RESREFC KE)=RKEINXRESREFC1) 

RESREFC EP)=EPIN*RESREFC1) 

I ERRP=L SWEEP/ 10 
KOUTPT  =  LSWEEP/1 0 


-  GROUP  40  SPECIAL  DATA  :  LOGICC 1 . . 10 ) , INTGRC 1 . . 10 ) , REC 21 . . 30 ) , 

NLSP<1>,NISP<1>,NRSP<1>, SPDATA< .F.>,LSPDAC1),ISPDAC1) ,RSPDAC1) 
USE  FIRST  10  ELEMENTS  OF  ARRAYS  LOGIC  8  INTGR  AND  21ST 
TO  30TH  OF  ARRAY  RE  FOR  TRANSFERRING  SPECIAL  DATA  FROM 
SATELLITE  TO  GROUND,  BUT  IF  REQUIREMENTS  EXCEED  THIS 
PROVISION  SET  SPDATA  =  .T.,  AND  DIMENSION  ARRAYS  LSPDA, 

ISPDA,  RSPDA  ABOVE  AND  IN  GROUND  AS  NEEDED,  AND  SET  HERE 


—  GROUP  42  RESTARTS  AND  DUMPS  1  SAVEM< . F . >, RESTRT< . F . > , KINPUT 
SAVEM= .TRUE . 

BFPLOT  = . TRUE. 

RESTRT  = . TRUE . 


-  GROUP  43  GRAFFIC  1 

GRAPHS <. F. >, 0RTH0G< . T . >, ANTSYM, NPRT<l>,ITITL<5x4H*x*x> 

-  FOR  A  GRAFFIC  RUN,  DIMENSION  PHI1  8  PHI2  AS  FOLLOWSt 

PHI1CNXXNYXNZXNM) 

PHI2CCNX+2)XCNY+2)XCNZ+2)XCNM+IBLK))  ,  WHERE 
NM=N0 .  OF  VARIABLES  STORED  +  DENSITYC -IES ) 

IBLK=0  IF  BLOCK= . FALSE. , =4  IF  A  3D  RUN, 

=3  IF  A  2D.YZ  RUN. 


IFCIRUN.EQ.l)  GO  TO  900 
900  CONTINUE 
C—  ALL  RUNS 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  USER  SECTION  3  ENDS. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  4  STARTS  1 


VAN05050 

VAN0506  0 

VAN0507  0 

VAN05080 

VAN05090 

VAN05100 

VAN05110 

VAN05120 

VAN05130 

VAN05140 

VAN05150 

VAN0516  0 

VAN05170 

VAN05180 

VAN05190 

VAN052Q0 

VAN05210 

VAN05220 

VAN05230 

VAN05240 

VAN05250 

VAN05260 

VAN0527  0 

VAN05280 

VAN05290 

VAN05300 

VAN05310 

VAN05320 

VAN05330 

VAN05340 

VAN05350 

VAN05360 

VAN0537  0 

VAN05380 

VAN05390 

VAN05400 

VAN05410 

VAN05420 

VAN05430 

VAN05440 

VAN05450 

VAN0546  0 

VAN0547  0 

VAN05480 

VAN05490 

VAN05500 

VAN0551 0 

VAN05520 

VAN05530 

VAN 05540 

VAN05550 

VAN0556  0 

VAN0557  0 

VAN05580 

VAN0559  0 

VAN05600 

VAN056 10 

VAN05620 

VAN05630 

VAN05640 

VAN05650 

VAN0566  0 

VAN0567  0 

VAN05680 

VAN05690 

VAN057  00 

VAN057 1 0 

VAN05720 

VAN05730 

VAN05740 

VAN05750 

VAN0576  0 
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FILEi  VAN4SAT  FORTRAN  A1 


WRITE  GENERAL  DATA  ON  TO  THE  GUSIE1.DTA  TAPE,  ETC... 

IFCSPDATA)  CALL  WRTSPC( LSPDA, NLSP, ISPDA, NISP, RSPDA, NRSP) 
IF(BLOCK)  CALL  WRTPORCPE, PN, PH, PC, NX, NY, NZ, IPLANE) 
IF(BFC)  CALL  WRTBFC(14,NBFC,XE,XW,YN,YS,ZH,ZL, 

&ND, NX+1 , NY+1 , NZ+1 , NZ, PRTBFC) 

OLD  PRACTICES  RETAINED  FOR  REFERENCE. 

IF(SPDATA)  CALL  SPCDAT(IRUN) 

IF( BLOCK)  CALL  PORDAT ( IRUN) 

IF(GRAPHS)  CALL  SORT(IRUN) 

IF(RESTRT)  GO  TO  902 
DO  901  INDVAR=1,25 

IF( IFIX( FIINITt INDVARHO .  1).NE. 10101)  GO  TO  901 
CALL  FLDDAT(IRUN) 

GO  TO  902 

901  CONTINUE 

902  CALL  DATAIOCWRT ,10) 

IFCMONITR)  CALL  DATAIO(WRT,-6) 

999  CONTINUE 
STOP 
END 

C*xx  IGEN=1  SO  BFCXYZ  NOT  REQUIRED. 

CXXX  COMMENT  OUT  BOTH  VERSIONS. 

C - 

SUBROUTINE  BFCXYZ  CNXP1,NYP1,NZP1) 

RETURN 

END 


VAN05770 
VAN05780 
VAN05790 
VAN05800 
VAN05810 
VAN05820 
VAN05830 
VAN05840 
VAN05850 
VAN05860 
VAN05870 
VAN05880 
VAN05890 
VAN05900 
VAN059 1 0 
VAN05920 
VAN05930 
VAN05940 
VAN05950 
VAN0596  0 
VAN05970 
VAN05980 
VAN05990 
VAN06000 
VAN06010 
VAN06020 
VAN06030 
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Ground  Listing 
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FILE*  VANTGRD  FORTRAN  Al 


c$directivexxmain  ami  leitner 

c  LECGRD  LAST  GEO.  NZ=27  NY=18  LAMINAR  FLOW 
C  XFILE  NAME:  MODBFCGD . FTN 

C  XINCLUDE  DED  SUBROUTINES:  THE  MODELS  OF  MAIN,  GROUND  8  STRIDE. 
C  XDOCUMENTATION:  PHOENICS  INSTRUCTION  MANUAL  (SPRING  1983) 

C  WITH  BODY-FITTED  COORDINATES  INSTRUCTION  SUPPLEMENT 

C  (SUMMER  1984). 

C  XSATELLITE  FILE  NAME:  MODSTL . FTN 

COMMON/ 1 SH I  FT/ 1 1 1 ( 57 ) , NFMAX 

C  SET  F-ARRAY  DIMENSION  AS  NEEDED,  &  SET  NFMAX  ACCORDINGLY. 

C  FOR  BFC'S  ALSO  SET  Fl-ARRAY  DIMENSION  AS  NEEDED  , AND  SET 
C  NF1MAX  ACCORDINGLY. 

COMMON/ FOB/ FI (10000) 

COMMON/ NF0B/NF1 MAX 
COMMON  F( 25000 ) 

NFMAX=25000 
NF1MAX=1 0000 
CALL  MAIN1 
STOP 
END 

CSDIRECTIVEXXGROUND 

SUBROUTINE  GROUND ( I RN, I CHAP, ISTP, ISWP, I ZED, INDVAR) 

INCLUDE  ( CMNGUS ) 

INCLUDE  ( GUSSEQ ) 

C  INCLUDE  NMLIST 

LOGICAL  BFC 

EQUIVALENCE  ( LOGIC( 20 ) , BFC  ) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  1  STARTS: 

C - 

C+++++MEANING  OF  SUBROUTINE  ARGUMENTS: 

C  IRN=RUN  NUMBER;  ICHAP=CHAPTER  CALLED;  ISTP=TIME  STEP; 

C  I$WP=SOLUTION  SWEEP;  IZED=Z-SLAB;  INDVAR:  SEE  CHAPTERS  BELOW. 
C+++++USER- INTRODUCED  VARIABLES  8  ARRAYS: 

C  TO  AVOID  CONFLICT  WITH  VARIABLE  NAMES  USED  IN  COMMON,  ALL 
C  VARIABLES  INTRODUCED  BY  THE  USER  SHOULD  HAVE  NAMES  STARTING 
C  WITH  *G'  IF  REAL,  'J'  IF  INTEGER,  AND  *Gf  OR  'J'  IF  LOGICAL. 

C  THUS  GDZ(IZ)  MIGHT  BE  A  Z-INTERVAL  ARRAY; 

C  GW1 ( IY, IX )  A  2-D  ARRAY  FOR  AXIAL  VELOCITY;  ETC. 

C  USER-GENERATED  SUBROUTINES  SHOULD  BE  NAMED  CORRESPONDINGLY, EG 

C  SUBROUTINE  GVISC(GTEMP, GCNC, GVSC) ,  FOR  COMPUTING  VISCOSITY 
C  FROM  CONCENTRATION  8  TEMPERATURE. 

C+++++GROUND-TO- EARTH  CONNECTING  SUBROUTINES: 

C  XUSE  GET ( NAME, G ARRAY, NY, NX)  TO  PUT  VALUES  OF  VARIABLE  NAMED 
C  1  NAME1  INTO  ARRAY  'GARRAY'  DIMENSIONED  GARRAY( NY, NX) . 

C  XUSE  SET(NAME,IXF,IXL,IYF,IYL, GARRAY, NY, NX)  TO  SET  VARIABLE 
C  1  NAME'  TO  GARRAY( IY, IX)  OVER  THE  REGION:  IXF-IXL  8  IYF-IYL . 

C  XUSE  PRN$LB( NAME)  TO  PRINT  VARIABLE  'NAME*  OVER  X-Y  PLANE. 

C  XUSE  ADD(NAME,IXF,IXL,  IYF,  IYL  ,  TYPE,  CM,  VM,  CVAR,  WAR,  NY,  NX) 

C  TO  ADD  SOURCE  TO  VARIABLE  NAMED  'NAME'  (SEE  CHAPTER  5). 

C  XUSE  READIZ(IZED)  IN  CHAPTERS  1,  2,  8,  8  9  TO  ACCESS  PI,. 

C  8  VOL , . . . AHDZ .  (SEE  FOOTNOTE  TO  LEGALITY  TABLE) 

C  XUSE  GET1D(NAME, GARRAY, NDIM)  TO  PUT  VARIABLE  NAMED  'NAME' 

C  ONE-D  ARRAY  ’GARRAY1  DIMENSIONED  NDIM,  THUS: 

C  CALL  GET 1 D( NAME, GNX, NX)  FOR  XG,...DXG  8  DIMENSION  GNX(NX); 

C  CALL  GET1D(NAME,GNY , NY)  FOR  YG,...RV  8  DIMENSION  GNY(NY) ; 

C  CALL  GET1D( NAME, GNZ, NZ)  FOR  ZG, . . . WGRID  8  DIMENSION  GNZ(NZ). 

C+++++LEGALITY  TABLE  FOR  USE  OF  EARTH-CONNECTING  SUBROUTINES: 


DM 


IN 


C 

C 

C 


ENTRIES  IN  TABLE  GIVE  CHAPTERS  IN  WHICH  SUBROUTINES  CAN  BE 
USED  FOR  VARIABLES  IN  LEFT-HAND  COLUMN.  (SUBROUTINE 
STRIDE  IS  REGARDED  AS  BEING  IN  CHAPTER  3) 


VAN0001 0 
VAN00020 
VAN00030 
VAN00040 
VAN00050 
VAN0006  0 
VAN00070 
VAN00080 
VAN00090 
VAN00100 
VAN00110 
VAN00120 
VAN00130 
VAN00140 
VAN00150 
VAN00160 
VAN00170 
VAN00180 
VAN00190 
VAN00200 
VAN0021 0 
VAN00220 
VAN00230 
VAN00240 
VAN00250 
VAN0026  0 
VAN00270 
VAN00280 
VAN00290 
VAN00300 
VAN00310 
VANO  0320 
VAN00330 
VAN00340 
VAN00350 
VAN0036  0 
VAN0037  0 
VAN00380 
VAN00390 
VAN00400 
VAN0041 0 
VAN  00420 
VAN00430 
VAN00440 
VAN00450 
VAN0046  0 
VAN0047  0 
VAN00480 
VAN00490 
VAN00500 
VAN0051 0 
VAN00520 
VAN00530 
VAN00540 
VANO  0550 
VAN0056  0 
VAN0057  0 
VAN00580 
VAN00590 
VAN006  00 


C 

c  • 

:  VARIABLE 
: 

: :  GET  8  : 

::  PRNSLB  : 

SET  : 

: 

ADD 

: 

: 

READIZ  : 

: 

GET1D  : 

: 

VMHUUOl U 

VAN00620 
VANO  06  30 
VAN00640 
VAN00650 

L 

C 

: P 1  -  RZ 

: :  ALL  : 

6  8  7: 

5 

: 

1,2, 8, 9: 

NONE  : 

c 

: P 1 0  -  RZH 

::3-7,  10-16: 

3  : 

NONE 

: 

NONE  : 

NONE  : 

VAN0066  0 

c 

: VOL  -AHDZ 

::  ALL  : 

3  : 

NONE 

: 

1,2, 8, 9: 

NONE  : 

VAN0067  0 

c 

:  D1  DP 

:  :  NONE  : 

10  : 

NONE 

: 

NONE  : 

NONE  : 

VAN00680 

c 

:D2DP 

:  :  NONE  : 

11  : 

NONE 

: 

NONE  : 

NONE  : 

VAN006  90 

c 

:  MU1 , MU1 H 

:  :  5,13-16  : 

12 

NONE 

: 

NONE  : 

NONE  : 

VAN007  00 

c 

:EXCO(L,H) 

:  :  NONE  : 

13  : 

NONE 

: 

NONE  : 

NONE  : 

VANO  07 10 

c 

:  CFP 

:  :  5  : 

14  : 

NONE 

: 

NONE  : 

NONE  : 

VAN00720 
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FILE* 

VANTGRD  FORTRAN  A1 

0 

i MDT  i i 

5  i 

■  15 

*  NONE  i 

i  NONE  i 

i  NONE 

0 

:  HST1 , HST2 1 i 

5  &  15  i 

>  16 

t  NONE  i 

i  NONE  i 

i  NONE 

C 

i XG  -WGRIDc 

NONE  i 

i  NONE 

NONE  i 

i  NONE 

i  ALL 

c  - 

c  NOTES  ON  ABOVE  TABLE* 

C  *IN  CHAPTERS  1,  2,  8,  8  9  VARIABLES  PI... DM  S  GEOMETRY 
C  VOL...AHDZ  CAN  BE  ACCESSED  BUT  ONLY  IN  CONJUNCTION  WITH 
C  USE  OF  READIZ,  THUS: 

C  DO  1  IZED=1 , NZ 

C  CALL  READIZC IZED) 

C  1  CALL  GET (  .  .  .  AS  REQUIRED..) 

C  ^GEOMETRY  ACCESSED  BY  READIZ  IS  THAT  AT  INITIAL  TIME. 

C  *D1DP  &  D2DP  ONLY  ACCESSIBLE  IN  UNSTEADY  FLOWS; 

C+++++GROUND  SERVICE  SUBROUTINES* 

C  *USE  CONTUR(NAME, IPLANE, ILOC, NINT, II, 12, Jl, J2, GARRAY, NDIM)  FOR 
C  LINE-PRINTER  PLOTS  OF  CONTOURS.  9  NAME 1  =  U1,...C4; 

C  9 IPLANE ' =  XPLANE,  YPLANE,  OR  ZPLANE;  ILOC  SETS  IX,  IY,  OR 

C  IZ  LOCATION  OF  IPLANE;  II,  12,  Jl,  S  J2  SET  FIRST  8  LAST 

C  CELLS  IN  HORIZ.  8  VERT.  ON  PLOT;  GARRAY  IS  1-D  WORKING  ARRAY 

C  OF  DIMENSION  NX*NY,  NX*NZ,  OR  NY*NZ  DICTATED  BY  IPLANE;  8 

C  NDIM  SETS  VALUE  OF  DIMENSION  OF  GARRAY. 

C  *USE  FLD2DAC TIT L E, GARRAY , NY, NX)  TO  PRINT  ANY  ARRAY  DIMENSIONED 
C  GARRAY ( NY, NX) ;  SET  'TITLE'  TO  REQUIRED  NAME  (  4  HOLLERITH 

C  CHARACTERS  ONLY) . 

C  *USE  FLD3DACTITLE, GARRAY, NX, NY, NZ, IPLANE, ILOC)  TO  PRINT  ANY 
C  ARRAY  DIMENSIONED  GARRAY (NX, NY, NZ)  IN  PLANE  SPECIFIED  BY 

C  'IPLANE'  8  'ILOC'  AS  FOR  CONTUR  ABOVE;  SET  'TITLE'  AS  FOR 

C  FLD2DA. 

C  VARIABLE  NAMES  FOR  USE  IN  GROUND* 

COMMON/TYPE/CELL, EAST, WEST, NORTH, SOUTH, HIGH, LOW, VOLUME, WALL 
COMMON/ VAR/ PI ,PP,U1,U2,V1,V2,W1,W2,R1,R2,RS, 
8KE,EP,H1,H2,H3,C1,C2,C3,C4,RX,RY,RZ,S1,S2 
COMMON/VAROL D/PIO , PP0,U10,U20,V10,V20,WI0, W20, RIO , R20 ,  RSO, 
8KE0,EP0,H10,H20,H30,C10,C20,C30,C40, RXO, RY0,RZ0,S10,S20 
COMMON/ VARL0W/PIL,PPL,U1L,U2L, V1L,V2L,W1L,W2L,R1L,R2L,RSL, 
8KEL,EPL,H1L,H2L,H3L,C1L,C2L,C3L,C4L,RXL,RYL,RZL,S1L,S2L 
C0MM0N/VARHI/P1H,PPH,U1H,U2H,V1H, V2H, W1H, W2H, R1H, R2H, RSH, 
8KEH,EPH,H1H,H2H,H3H,C1H,C2H,C3H,C4H, RXH, RYH, RZH , S1H , S2H 
COMMON/GMTRY/VOL, VOL 0 , AEAST , ANORTH , AHIGH,AEDX, ANDY,  AHDZ 
COMMON/PROP/ D1,D2,D1 DP, D2DP , MUI , MUI L AM, EXCO ,CFP , MDT, HST1 , HST2 
COMMON/PRPOL D/ DIO , D20 
COMMON/ P RP L  OW/ D1 L , D2L , EXCOL 
COMMON/ PRPH 1/ D1H , D2H, MU1H, EXCOH 
COMMON/VARNX/XG, XU , DXU, DXG 
COMMON/VARNY/YG, YV , DYV, DYG, R, RV 
COMMON/ VARNZ/ZG, ZW1 , DZW, DZG, WGRID 
COMMON/ GDMSCI/XP LANE ,YP LANE, ZPLANE, ITNO 
COMMON/GDMSCL/LSLAB, MSLAB, HSLAB, L AMMU 
REAL  NORTH, LOW 

INTEGER  P1,PP,U1,U2,V1,V2,W1,W2,R1,R2,RS, 

8EP,H1,H2,H3,C1,C2,C3,C4,RX,RY,RZ,S1,S2 
INTEGER  P10,PP0,U10,U20,VI0,V20,W10, W20 ,R10,R20,RS0, 

8EP0,H10,H20,H30,C10,C20,C30,C40,RX0,RY0, RZO, SIO, S20 
INTEGER  PIL,PPL,U1L,U2L,V1L, V2L , W1 L , W2L , RI L , R2L , RSL , 

8EPL,H1L,H2L,H3L,C1L,C2L,C3L,C4L, RXL ,RYL,RZL,S1L,S2L 
INTEGER  PIH,PPH,U1H,U2H, V1H, V2H, WIH, W2H, RIH, R2H, RSH, 

8EPH,H1H,H2H,H3H,C1H,C2H,C3H,C4H, RXH, RYH, RZH,S1H,S2H 
INTEGER  VOL, VOL 0, AEAST, ANORTH, AHIGH,AEDX, ANDY, AHDZ 
INTEGER  D1 , D1DP , D2, D2DP , EXCO, CFP , HST1 , HST2 
INTEGER  D10,D20,D1L,D2L, EXCOL , D1H, D2H, EXCOH 
INTEGER  XG, XU, DXU, DXG,YG,YV, DYV, DYG, R, RV, ZG,ZW1 , DZW, 

8DZG, WGRID 

INTEGER  XPLANE,YPLANE, ZPLANE 
LOGICAL  LSLAB, MSLAB, HSLAB, LAMMU,LSPDA 
EQUIVALENCE  (Ml , Rl ) , ( M2, R2) 

C  SATLIT-EQUIVALENT  IRUN* 

EQUIVALENCE  ( IRUN , INTGR( 11 ) ) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  1  ENDS. 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  USER  SECTION  1  STARTS* 

C  ARRAYS  (  DIMENSIONED  NY, NX  )  FOR  USE  WITH  'ADD'* 

DIMENSION  CV AR(1,1),VVAR(1,1), CM (1,1),VM(1,1),ZER 0(1,1) 
DIMENSION  GP(30,1),GH(30,1),GD(30,1),GV(30,1),GW(30,1) 
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VAN007  30 
VAN00740 
VAN00750 
VAN0076  0 
VAN0077  0 
VAN00780 
VAN00790 
VAN00800 
VAN0081 0 
VAN00820 
VAN00830 
VAN00840 
VAN00850 
VAN0086  0 
VAN00870 
VAN0088  0 
VAN0089  0 
VAN009  00 
VAN00910 
VAN00920 
VAN00930 
VAN00940 
VAN00950 
VAN0096  0 
VAN0097  0 
VAN00980 
VAN00990 
VAN01000 
VAN01010 
VAN01020 
VANO 1030 
VANO 1040 
VAN01050 
VANO 1 06  0 
VANO 1 07  0 
VANO 1080 
VAN01090 
VANO 1 100 
VAN01110 
VAN01120 
VAN01130 
VANO 1 140 
VAN01150 
VANO 1 1 6  0 
VANO 1 170 
VAN01180 
VANO 1 190 
VANO 1200 
VANO 1210 
VAN01220 
VAN01230 
VAN 01240 
VAN01250 
VANO 126  0 
VAN01270 
VANO 1280 
VANO 1 29  0 
VANO 1300 
VAN01310 
VAN01320 
VANO 1 330 
VAN01340 
VAN01350 
VAN0136  0 
VAN0137  0 
VAN01380 
VAN01390 
VAN01400 
VANO 14 1 0 
VAN01420 
VAN01430 
VAN01440 
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1  ,GMACH(30,1),GTEMP(30,1),GVISC(30,1),GWH(30,1),GWM(30,1)  VAN01450 

2  ,GKE(30,1),GC3(30,1),GYG(30,1),GXX(30,1),GYY(30,1),GZZ(30,1)  VAN01460 

C  SPECIAL-DATA  ARRAYS  DIMENSIONED  8  DIMENSION  VALUES  SET  HERE:  VAN01470 

DIMENSION  LSPDA(1),ISPDA(1),RSPDA(1)  VAN01480 

C  USER  PLACES  HIS  VARIABLES,  ARRAYS,  EQUIVALENCES  ETC.  HERE.  VAN01490 

EQUIVALENCE  ( RAIR, RE(21 ) ) , (GAMA, RE(22) ) , (GSWP, RE( 23) ) ,  VAN01500 

1CGPR, REC  24)),(GTW,RE(25)), (GEMU1 , RE( 26 ) ) , ( JEMU1 ,INTGR(1))  VAN01510 

DATA  NLSP, NISP, NRSP/1 , 1 , 1/  VAN01520 

DATA  CVAR,  WAR,  CM,  VM,  ZERO/5XO  .  0/  VAN01530 

C  USER  PLACES  HIS  DATA  STATEMENTS  HERE.  VAN01540 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  USER  SECTION  1  ENDS.  VAN01550 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  2  STARTS:  VAN01560 

C  PLEASE  DO  NOT  ALTER,  OR  RE-SET,  ANY  OF  THE  REMAINING  VAN01570 

C  STATEMENTS  OF  THIS  SECTION.  VAN01580 

DATA  NUMCH4  /  0  /  VAN01590 

IFCSPDATA)  VAN01600 

8CALL  RDSPCC IRN,INTGR(12),LSPDA,NLSP,ISPDA,NISP,RSPDA,NRSP)  VAN01610 

CALL  GRDUTYC IRN, ICHAP, I ZED, INDVAR)  VAN01620 

IFCBFC)  CALL  BFCGRDC IRN, ICHAP, ISWP, IZED, INDVAR)  VAN01630 

IFC ICHAP . EQ . -5)  GO  TO  10  VAN01640 

IF( ICHAP . L E . 0 . OR . ICHAP . GT . 16 )  RETURN  VAN01650 

GO  TO  (100,200,300,4999,500,600,700,800,900,1000,1100,1200,  VAN01660 

81300, 1400, 1500, 1600), ICHAP  VAN01670 

RETURN  VAN01680 

4999  NUMCH4=  NUMCH4  +  1  VAN01690 

IF  (MOD( NUMCH4, 2) . EQ . 1 )  GO  TO  400  VAN01700 

RETURN  VAN01710 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  STANDARD  SECTION  2  ENDS.  VAN01720 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  USER  SECTION  2  STARTS:  VAN01730 

c -  VAN01740 

C  CHAPTER  0:  MODIFY  SATLIT  DATA,  AT  START  OF  EACH  IRN.  VAN01750 

c -  VAN01760 

10  CONTINUE  VAN01770 

C  IF( .NOT.NAMLST)  RETURN  VAN01780 

C  IF(IRN.EQ.NRUN)  DATFIL= . FALSE.  VAN01790 

C—  READ  SATLIT  DATA  NAMELIST  HERE  VAN01800 

C  CALL  WRIT40(40HENTER  NAMELIST  DATA  FOR  GROUPS  1  TO  24  )  VAN01810 

C  READ( 20 , G1G24 )  VAN01820 

C  CALL  WRIT40(40HENTER  NAMELIST  DATA  FOR  GROUPS  25  TO  42  )  VAN01830 

C  READ(20,G25G42)  VAN01840 

RETURN  VAN01850 

c -  VAN01860 

C  CHAPTER  1:  CALLED  AT  THE  START  OF  EACH  TIME  STEP.  VAN01870 

C  SET  'DT'  HERE  WHEN  TLAST  SET  NEGATIVE  IN  BLOCK  DATA.  VAN01880 

C  ' ATIME  +  DT*  GIVES  THE  END  TIME  OF  THE  CURRENT  TIME  STEP.  VAN01890 

C  NOT  ACCESSED  IF  STEADY, OR  PARABOLIC.  VAN01900 

C -  VAN01910 

100  CONTINUE  VAN01920 

RETURN  VAN01930 

C -  VAN01940 

C  CHAPTER  2:  CALLED  AT  THE  START  OF  EACH  SWEEP.  VAN01950 

c -  VAN01960 

200  CONTINUE  VAN01970 

RETURN  VAN01980 

C -  VAN01990 

C  CHAPTER  3:  CALLED  AT  THE  START  OF  EACH  SLAB;  VAN02000 

C  NOT  ACCESSED  IF  PARABOLIC,'  BUT  'STRIDE'  IS.  VAN02010 

C -  VAN02020 

300  CONTINUE  VAN02030 

RETURN  VAN02040 

c -  VAN02050 

C  CHAPTER  4:  CALLED  AT  THE  START  OF  EACH  RE-CALCULATION  OF  VAN02060 

C  VARIABLES  P1,...C4  AT  CURRENT  SLAB.  ITNO=  ITERATION  NUMBER.  VAN02070 

C — ; -  VAN02080 

400  CONTINUE  VAN02090 

RETURN  VAN02100 

c -  VAN02110 

C  CHAPTER  5:  GROUND  CALLED  WHEN  SOURCE  TERM  IS  COMPUTED.  VAN02120 

C  INDVAR  GIVES  DEPENDENT  VARIABLE  IN  QUESTION  IE.  U1,...C4.  VAN02130 

C  TO  ADD  SOURCE  TO  DEPENDENT  VARIABLE  Cl(SAY)  FOR  IX=IXF,IXL  VAN02140 

C  AND  IY=IYF , IYL  INSERT  STATEMENT:  VAN02150 

C  IF(INDVAR.EQ.Cl)  VAN02160 
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FILE*  VANTGRD  FORTRAN  A1 


&CALL  ADD(INDVAR,IXF,IXL,IYF,IYL,TYPE,CM,VM,CVAR,VVAR,NY,NX) 

NOTES  ON  'ADD1 i 

XS0URCE=  (CVAR(IY, IX)+AMAX1(0 . 0, MASFLO) )x(VVAR( IY, IX)-PHI ) , 

WHERE  'PHI*  IS  IN-CELL  VALUE  OF  VARIABLE  IN  QUESTION, 
x 'MASFLO ' 3  CM(IY,IX)x(VM(IY,IX)-P), 

WHERE  'P'  IS  THE  IN-CELL  PRESSURE. 

XFOR  INDVAR=  Ml,  OR  =M2,  SOURCE  ADDED  IS  'MASFLO '  ONLY, 

EXCEPT  FOR  0NEPHS=.F.  8  MASFLO  <  0.0  (IE.  OUTFLOW)  WHEN 
CMC IY, IX)  IS  MULTIPLIED  BY  R1*D1  (FOR  Ml)  a  R2XD2  (FOR  M2). 

XBOTH  'CVAR 1  &  'CM*  ARE  MUTLIPLIED  BY  CELL-GEOMETRY  QUANTITY 
DICTATED  BY  SETTING  OF  ’TYPE 1  (=CELL,  EAST  AREA, .. VOLUME ) . 
XTYPE-SPECIFIED  AREAS  ARE  CALCULATED  AS  IF  BLOCKAGE  ABSENT, 

BUT  ’VOLUME’  WITH  ACCOUNT  FOR  ITS  PRESENCE. 

XFOR  ALL  SOLVED  VARIABLES,  INCLUDE  DING  Ml  (  a  M2  WHEN  QNEPHS=F), 
IF  ' CM ' >  0.0  CALL  ’ADD’;  FOR  Ml  a  M2  ALTHOUGH  ’CVAR’  a  ’WAR’ 
HAVE  NO  SIGNIFICANCE  THEY  MUST  BE  ENTERED  AS  ARGUMENTS. 

X 'CVAR 1 ,  'WAR',  'CM'  a  'VM'  MUST  BE  DIMENSIONED  NY, NX. 

500  CONTINUE 
RETURN 

CHAPTER  6 i  CALLED  AT  THE  END  OF  EACH  VARIABLE-RECALCULATION 
CYCLE  COMMENCED  AT  CHAPTER  4.  ITNO  =  ITERATION  NUMBER. 

600  CONTINUE 
RETURN 

CHAPTER  7 i  CALLED  AT  END  OF  EACH  SLAB-WISE  CALCULATION. 

700  CONTINUE 

IF( FLOAT ( ISWP ) . LT . GSWP)  RETURN 
CALL  GET( PI , GP, NY, NX) 

CALL  GET( HI , GH, NY, NX) 

CALL  GET( D1 ,GD, NY, NX) 

CALL  GET( VI , GV, NY, NX) 

CALL  GET (W1 , GW, NY, NX) 

CALL  GET(KE,GKE, NY, NX) 

C  CALL  GET1D(YG,GYG,NY) 

CALL  GREDK39,  IZED,  GYG,  NY,  NX) 

CALL  GRED3( 57 , IZED, GXX, GYY , GZZ, NY, NX ) 

GCP=RAIR/(1. -1/GAMA) 

DO  701  1=1, NY 

GSON=SQRT(GAMAXGP(I, 1)/GD(I,1)) 

GAV=SQRT(GV( I , 1 )XX2+GW( I, l)xx2) 

GMACH( 1,1) =GAV/GSON 

701  GTEMP(I, 1)=GP(I,1)/GD(I,1)/RAIR 

C  701  GTEMP(I>l)=(GH(I,l)-GW(I,l)xx2/2.-GV(I,l)XX2/2.)/GCP 
CALL  SET ( Cl , 1 , NX , 1 , NY, GMACH, NY, NX ) 

CALL  SET(C2,1,NX,1,NY, GTEMP, NY, NX) 

C - CALCULATE  DY1  CF  ST  H(CONVECTIVE  COEF.)  Q  TAU  TR 

I F( J EMU1 . NE . 2 )  GOTO  702 

C - TURBULENT  VALUES 

GCF=2./GW(NY, 1)XX2XGKE(1, l)/3.33xGD(l, 1)/GD(NY, 1) 

C7  GCF=GCFXGD(NY, 1)/GD(1, 1)XGTEMP(NY, 1)/GTEMP(1, 1)XGP(1, 1)/GP(NY, 1) 
GST=GCF/2./GPRXX.666 
GHH=GD(NY, 1)XGCPXGW(NY, 1)XGST 
GR=GPRXX . 333 

GTR=GTEMP( NY, 1 ) X( 1 . +GRx(GAMA-l . )/2 . xGMACH(NY, 1 )xx2) 

C  1(1 .+( GAMA-1 . )/2.xGMACH(NY, 1)XX2) 

GQ=GHHX(GTR-GTW) 

GOTO  703 

C - LAMINAR  VALUES 

702  CONTINUE 

IF( JEMU1 .EQ.-l)  GEMU1=GVISC( 1,1) 
GQ=GEMU1/GPRX(GH(1,1)-GTWXGCP)/GYG(1,1) 

GR=GPRxx .5 

GTR=GTEMP(NY,1)X(1 .+GRX(GAMA-1. )/2.XGMACH(NY,l)xx2) 

C  1(1 .+(GAMA-1 . )/2.xGMACH(NY, 1)XX2) 

GHH=GQ/(GTR-GTW) 

GST  =  GHH/ ( GD ( NY , 1 ) XGW ( NY , 1 ) XGCP ) 

GTAU=GEMUlxGW( 1 , 1 )/GYG( 1,1) 

GCF=GTAUX2./(GD(NY,1)XGW(NY,1)XX2) 


VAN02170 
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FILEi  VANTGRD  FORTRAN  A1 


703  GC3( 1 ,  1 )=GYG( 1 , 1 ) 

GC3(2,1)=GCF 
GC3(3,1)=GST 
GC3(4, l)=GCF/2./GST 
GC3(5,1)=GHH 
GC3(6,1)=GQ 
GC3(7,1)=GTAU 
GC3(8, 1  )=GTR 
GC3( 9 , 1 ) =GTR-GTW 

GC3(10,1)=GD(NY,1)*GW(NY,1)XGZZ(1,1)/GEMU1 
GC3( 11 , 1 )=GZZ( 1, 1 ) 

GC3(12,1)=GEMU1 

GC3(13,l)=GD(NY,l)*GW(NY,l)XGYG(l,l)/GEMUlxSQRT(ABS(GCF/2.)) 
CALL  SET(C3,1,NX,1,NY,GC3,NY,NX) 

RETURN 


CHAPTER  8:  CALLED  AT  THE  END  OF 
NOT  ACCESSED  IF  PARABOLIC. 

EACH 

SWEEP; 

800 

CONTINUE 

RETURN 

CHAPTER  9:  CALLED  AT  THE  END  OF 
NOT  ACCESSED  IF  PARABOLIC. 

EACH 

TIME  STEP; 

900 

CONTINUE 

RETURN 

CH^TER  10>  SET  PHASE  1  DENSITY 
SET  CURRENT-Z  'SLAB'  DENSITY,  D1 
EG.  IF(MSLAB)  CALL  SET(D1,1,NX, 

HERE 
,  IF 

1 ,  NY, 

WHEN  IRH01=-1  IN  DATA. 
MSLAB= .T . , 

GDI, NY, NX) . 

SET  NEXT  LARGER-Z  'SLAB'  DENSITY,  D1H,  IF  HSLAB=.T.  &  PARAB=F 
EG.  IF(HSLAB)  CALL  SETCD1H, 1 , NX, 1 , NY, GD1H, NY, NX) . 

SET  D(LN(D1))/DP  (IE.  D1DP)  FOR  UNSTEADY  FLOW, 

EG.  IF(MSLAB)  CALL  SET ( D1 DP, 1 , NX, 1 , NY, GD1DP, NY, NX) . 


1000  CONTINUE 

IF  (MSLAB)  GO  TO  101 

JP1=P1H 

JH1=H1H 

JD1=D1H 

JW1=W1H 

JV1=V1H 

GO  TO  102 

101  JP1=P1 
JH1=H1 
JD1=D1 
JW1=W1 
JV1=V1 

102  CALL  GET ( JP1 , GP, NY, NX) 

CALL  GET( JH1 , GH, NY, NX) 

CALL  GET( JW1 ,GW, NY, NX) 

CALL  GET( JV1 , GV , NY, NX) 

IF(IZED.EQ.l)  GOTO  105 
IF(IZED.EQ.NZ)  GOTO  109 

C - IZED=2,  NZ-1 

DO  103  IX=1 , NX 

DO  103  IY=1 , NY 

IF(HSLAB)  GOTO  104 

GWA=(GW(IY,IX)+GWMCIY,IX))/2. 

GWM(IY,IX)=GW(IY,IX) 

GOTO  115 

104  GWA=(GW(IY,IX)+GWH(IY,IX))/2. 

GWH( IY, IX)=GW( IY, IX) 

115  GHS=GH(IY,IX)-(GWAXX2+GV(IY,IX)X*2)/2. 
IF(GHS.LE.1.E5)  GHS=1.E5 

103  GD(IY,IX)=  GP( IY, IX )/( 1-1/GAMA )/GHS 
GOTO  113 

c - IZED=1 

105  DO  106  IX=1 , NX 
DO  106  IY=1 , NY 

GHS=GH( IY, IX)-(GW(IY,IX)XX2+GV(IY,  IX)**2)/2. 
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C 


C 


c 


107 

106 


109 


111 


112 

110 


IFCGHS . LE. 1 . E5)  GHS=1.E5 
GDCIY,IX)=  GPCIY, IX)/ C 1-1/GAMA )/GHS 
IFCHSLAB)  GOTO  107 
GWMCIY,IX)=GWCIY,IX) 

GOTO  106 

GWHCIY,IX)=GWCIY,IX) 

CONTINUE 
GOTO  113 

- IZED=NZ 


DO  110  IX=1,NX 
DO  110  IY=1 , NY 
IFCHSLAB)  GOTO  111 

GHS=GHCIY, IX)-CGWMC IY , IX) XX2+GVC IY, IX)xx2)/2.  * 
IFCGHS. LE.1.E5)  GHS=1.E5 
GWMCIY,IX)=GWCIY,IX) 

GOTO  112 

GHS=GHCIY, IX)-CGWHC IY, IX )x#2+GVC IY, IX)xx2)/2 . 
IFCGHS. LE.1.E5)  GHS=1.E5 
GWHCIY, IX) =GWC IY , IX) 

GDC IY, IX )=  GPCIY, IX)/ Cl -1/GAMA )/GHS 
CONTINUE 


113  CONTINUE 

CALL  SETC J D1 , 1 , NX, 1 , NY, GD, NY, NX) 
RETURN 


CHAPTER  Hi  SET  PHASE  2  DENSITY  HERE  WHEN  IRH02=-1  IN  DATA. 
SET  CURRENT-Z  'SLAB'  DENSITY,  D2,  IF  MSLAB= . T . , 

EG.  IFCMSLAB )  CALL  SETC D2, 1 , NX, 1 , NY,GD2, NY, NX) . 

SET  NEXT  LARGER-Z  1  SLAB 1  DENSITY,  D2H ,  IF  HSLAB= . T .  a  PARAB  =  F 
EG.  IFCHSLAB)  CALL  SETC D2H, 1 , NX, 1 , NY,GD2H, NY, NX) . 

SET  DC  LNC  D2 ) )/DP  FOR  UNSTEADY  FLOW, 

EG.  IFCMSLAB)  CALL  SETC D2DP, 1 , NX, 1 , NY, GD2DP, NY, NX) . 


1100  CONTINUE 
RETURN 


CHAPTER  1 2 i  SET  PHASE  1  VISCOSITY  HERE  WHEN  IEMU1=-1  IN  DATA. 
SET  CURRENT-Z  1  SLAB 1  VISCOSITY  CMU1),  IF  MSLAB= . T . , 

EG.  IFCMSLAB)  CALL  SETC MU1 , 1 , NX, 1 , NY, GVISC, NY, NX) . 

SET  NEXT  LARGER-Z  'SLAB*  VISC.  CMU1H),  IF  HSLAB= . T .  a  PARAB=F 
EG.  IFCHSLAB)  CALL  SETCMU1H, 1 , NX, 1 , NY, GVSCH, NY, NX) . 

CHAPTER  ALSO  ACCESSED  WHEN  EMULAM=-1.0  IN  DATA,  SO  THAT  THE 
LAMINAR  VISCOSITY  WHICH  APPEARS  IN  WALL  FUNCTIONS  a  IN  THE 
KE-EP  TURBULENCE  MODEL  CIEMU1=2)  MAY  BE  SET  NON-CONSTANT. 

SET  CURRENT-Z  'SLAB*  VALUE  CMU1LAM)  WHEN  LAMMU= . T . , 

EG.  IFCLAMMU)  CALL  SETCMU1LAM, 1 , NX, 1 , NY, GVSCL , NY, NX) . 


1200  CONTINUE 

GCP=RAIR/C1. -1/GAMA) 

IF  CHSLAB)  GOTO  122 
CALL  GETCH1,GH,NY,NX) 

CALL  GETC W1 , GW, NY, NX) 

CALL  GETC VI , GV, NY, NX) 

GOTO  123 

122  CALL  GET C  H1H, GH, NY, NX) 

CALL  GETC W1H, GW, NY,  NX) 

CALL  GETC V1H, GV, NY, NX) 

123  CONTINUE 

DO  121  IX=1 , NX 
DO  121  IY=1 , NY 

GTMP=CGHC IY, IX)-GWC IY, IX) XX2/2 . -GVC IY, IX )xx2/2 . )/GCP 
IFCGTMP.LT. 150.)  GTMP=150. 

121  GVISCCIY, IX) =1 .716E-05XCGTMP/273. )^^0 .666 
121  IFC  GVISCC IY, IX) .LE. .8E-5)  GVISCC IY, IX) = . 8E-5 
IF  CMSLAB)  CALL  SETC MU1 , 1 , NX, 1 , NY, GVISC, NY, NX) 

IF  CHSLAB)  CALL  SETC MU1H, 1 , NX, 1 , NY, GVISC, NY, NX) 

IF  CLAMMU)  CALL  SETC MU1LAM, 1 , NX, 1 , NY, GVISC, NY, NX) 
RETURN 


CHAPTER  13*  SET  EXCHANGE  COEFFICIENT  CE.C.)  FOR  VARIABLE 
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FILEi  vantgrd  fortran  ai 


INDVAR  WHEN  SIGMA(INDVAR)=-1 . 0  IN  DATA. 

SET  CURRENT-Z  'SLAB'  E.C.  (EXCO)  IF  MSLAB=.T., 

EG.  IF(MSLAB)  CALL  SETC EXCO, 1 , NX, 1 , NY, GEXCO, NY, NX) . 

SET  NEXT  SMALLER-Z  'SLAB'  E.C.  (EXCOL)  IF  LSLAB= . T . , 

EG.  IF(LSLAB)  CALL  SET( EXCOL , 1 , NX, 1 , NY, GEXCOL , NY, NX) . 

SET  NEXT  LARGER-Z  'SLAB'  E.C.  (EXCOH)  IF  HSLAB=.T., 

EG.  IF(HSLAB)  CALL  SET ( EXCOH, 1 , NX, 1 , NY , G EXCOH, NY ,  NX) . 

NOTE i  FOR  MSLAB,  INDVAR=U1, . .C4;  FOR  LSLAB,  INDVAR=U1L, . .C4L 
8  FOR  HSLAB,  INDVAR=U1H, . . C4H .  IF  PARAB= . T .  SET  MSLAB  ONLY. 


1300 

CONTINUE 

RETURN 

CHAPTER  14 i  SET 
WHEN  ICFIP  =  -1 
SPEED  OF  PHASES) 

INTER-PHASE  FRICTION  COEFFICIENT  CCFP)  HERE 

IN  DATA;  ITS  UNITS  =  FORCE  /  (CELL  X  RELATIVE 

1400 

CONTINUE 

RETURN 

CHAPTER  15:  SET 
HERE  WHEN  IMDOT 

INTER-PHASE  MASS-TRANSFER  RATE  PER  CELL  (MDT) 
=  -1  IN  DATA. 

1500 

CONTINUE 

RETURN 

CHAPTER  16 i  SET  HERE  PHASE  182  SATURATION  ENTHALPIES 
(  HST1  8  HST2)  WHEN  IHSAT  =  -1  IN  DATA. 

1600 

CONTINUE 

RETURN 

END 
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